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We study collective excitations of rotational and spin states of an ensemble of polar molecules, 
which are prepared in a dipolar crystalline phase, as a candidate for a high fidelity quantum memory. 
While dipolar crystals are formed in the high density limit of cold clouds of polar molecules under ID 
and 2D trapping conditions, the crystalline structure protects the molecular qubits from detrimental 
effects of short range collisions. We calculate the lifetime of the quantum memory by identifying the 
dominant decoherence mechanisms, and estimate their effects on gate operations, when a molecular 
ensemble qubit is transferred to a superconducting strip line cavity (circuit QED). In the case rota- 
tional excitations coupled by dipole-dipole interactions we identify phonons as the main limitation 
of the life time of qubits. We study specific setups and conditions, where the coupling to the phonon 
modes is minimized. Detailed results are presented for a ID dipolar chain. 



PACS numbers: 03.67.Lx, 33.80.Ps, 85.25.Cp, 61.50.-f 



I. INTRODUCTION 



In recent publications [l|, Q we have studied trapped 
polar molecules strongly coupled to a superconducting 
microwave strip line cavity, which represents a basic 
building block for hybrid quantum circuits, interfacing 
high-fidelity molecular quantum memory with solid state 
elements such as Cooper pair boxes (CPB) d, R IEJILj 
superconducting flux qubits J7|] or quantum dots |8l.l9l.lld. 

This suggests a hybrid quantum computing scenar- 
ios with the goal of combining the advantages of quantum 
optical and solid state implementations by interfacing 
molecular and solid state qubits in compatible experi- 
mental setups Polar molecules provide two key fea- 
tures for these interfaces. First, the long coherence times 
for qubits stored in polar molecules is based on iden- 
tifying long-lived rotational or electron and nuclear spin 
states in the electronic and vibrational ground state man- 
ifold. Second, the strong coupling of molecular qubits to 
the microwave cavity is based on transitions between ro- 
tational excitations (in the few GHz domain), with large 
electric dipole moments of the order of a few Debye. In 
this context, Ref. pj has studied the storage of single po- 
lar molecules on electric molecular chips ~ 100 nm above 
a superconducting strip line cavity, resulting in a single 
molecule - single photon coupling g of up to 1 MHz. Here, 
the strong coupling to the microwave cavity provides the 
additional features of cooling of the center-of-mass mo- 
tion of molecules and read out of the molecular qubits. 
In Ref. ^ a setup was studied (see Fig. [J), where a qubit 
was stored in the collective spin (or rotational) excitation 
of a thermal ensemble of N cold polar molecules, with the 
advantage of an enhanced collective coupling gjq = 
of the qubit to the cavity (typically of the order of 1 to 10 
MHz for TV = 10* — 10^ and trapping distances ^ 10/im), 
but at the expense of introducing dephasing of the qubit 
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Figure 1: a) Schematic picture of a hybrid quantum com- 
puter consisting of Cooper Pair Boxes (CPB) representing a 
solid state quantum processes and molecular ensemble qubits 
(MEQ) acting as long-lived quantum memories. Quantum in- 
formation is shuttled between to two systems via a supercon- 
ducting strip line cavity, b) Dipolar crystal of polar molecules. 
Under ID (or 2D) trapping conditions and dipole moments 
aligned by a strong electric field, repulsive dipole-dipole inter- 
actions stabilize the molecules against short-range collisions. 
See text for more details. 



due to state dependent collisions. 

In the present work we will investigate a molecular en- 
semble quantum memory in the form of a dipolar (self- 
assembled) crystal. Formation of a molecular dipolar 
crystal (MDC) is based on inducing electric dipole mo- 
ments fiind of polar molecules with an external DC field, 
which gives a I'-ind/'^'^^f^'''^ = C^/r^ interaction between 
molecules at distance r. For molecules confined to ID 
or 2D by an external trapping potential and dipole mo- 
ments aligned perpendicular, this interaction is repul- 
sive and allows the formation of a high-density crystal, 
where molecules perform small oscillations around their 
equilibrium positions, reminiscent of a Wigner crystal 
of trapped ions [31 . In this crystalline phase close en- 
counter collisions are strongly suppressed in comparison 
with thermal ensembles. 

The paper is organized as follows. We start in Sec. |TT] 
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with a brief review on molecular ensemble qubits and the 
hybrid quantum processor proposed in Ref. [4]. In Sec. Hill 
we then study the dynamics of collective rotational exci- 
tations in a MDC due to state dependent dipole-dipole 
interactions and the coupling to motional degrees of free- 
dom. We derive a general model for this system and 
calculate the resulting limitations for the lifetime of en- 
semble qubits stored in rotational degrees of freedom. In 
Sec.|lV]we extend our model to molecules with an addi- 
tional spin degree of freedom and show how we can realize 
a highly protected quantum memory using spin ensemble 
qubits in a MDC. In Sec. |V]we discuss in some detail a 
potential experimental implementation of a MDC using 
electrostatic traps, where in particular we study the ef- 
fects of an additional longitudinal confining potential. A 
summary and concluding remarks are given in Sec. IVII 

II. OVERVIEW AND BACKGROUND 
MATERIAL 

Before introducing our models of ensemble qubits in 
dipolar crystals in detail we find it useful to summarize 
briefly ensemble qubits and their coupling to microwave 
cavities, and dipolar crystals, to introduce and motivate 
the models for rotational and spin qubits described in the 
following sections. 

A. Molecular ensembles coupled to a 
superconducting strip line cavity 

We consider the setup of Fig. [T] a) , where a molecu- 
lar ensemble is coupled to a superconducting microwave 
cavity. In addition, the cavity could be coupled strongly 
to a CPB (or a quantum dot) representing a circuit QED 
system Q. As discussed in detail in Sec. IIII B[ molecular 
spectroscopy allows us to identify long-lived rotational 
states \g) and |e) within the electronic and vibrational 
ground state manifold. We assume that the ensemble 
of N polar molecules is prepared initially in the state 
I51 ■ ■ -ffAf) and is coupled to a single mode of a super- 
conducting microwave cavity with a frequency lOc close 
to the rotational transition frequency Ueg- The dipole 
coupling between the molecules and the cavity can then 
be written in the form 

-ffcav-mol = ngVN{Rlc + Re) , (1) 

with g the single molecule vacuum Rabi frequency and 
c (c^) the cavity annihilation (creation) operator. The 
collective molecular operator R\ creates symmetric Dicke 
excitations \ne) with the lowest two states 

|0e) = |5i...5jv), 

|le) = i?l|Oe) = l/V]V^|gi...e,...5jv), 

i 

representing an ensemble qubit, in addition to higher ex- 
citations of the form |2e) = l/\/2(i?t)2|0e), etc. For low 



number of rotational excitations the operator i?e fulfills 
approximate bosonic commutation relations [i?ei^e] — 1- 
As noted above, for typical experimental parameters the 
collectively enhanced coupling strength gjq = g\/N can 
be of the order of ten MHz, exceeding experimentally 
demonstrated decay rates of high-Q superconducting 
strip line cavities by several orders of magnitude [3, ■ 
For a single molecular ensemble coupled to a single 
CPB the total Hamiltonian for the hybrid system is 

Hsys = i?CPB + hUcC^C + hlOegRlRe 

0. J. a. V'^ / 

+ hgc{(J+c + o-_cT) + hgN{R\c + c^ Re) ■ 

Here terms in the first line represent the bare Hamilto- 
nian operators for the CPB, the cavity and the molec- 
ular ensemble. The second line of Eq. ^ contains the 
Jaynes-Cummings type interaction between cavity and 
CPB with a vacuum Rabi frequency gc (see Ref. jOl) as 
well as the molecule-cavity interaction given in Eq. ([1]). 
Interactions can in principle be controlled by tuning fre- 
quencies in and out of resonance or, as we explain later 
in this paper (see Sec. IVD[) . by using a switchable Ra- 
man process. Hamiltonian ^ is readily generalized to 
multiple molecular ensembles and CPBs. 

In summary Hamiltonian ^ provides the basic ingre- 
dient to convert the state of the CPB into a cavity pho- 
ton superposition state \ipc) = ct\Oc) +/3|lc) [ll], and in 
a second step to map it on an ensemble superposition 

(a|Oe) + /3|lc))|0e) ^ |0e)(a|0e) -|- /3|le)) , (3) 

and vice versa. As discussed in Ref. 2] coupling of a 
molecular ensemble qubit to higher excitations |2e) etc. 
can be suppressed by employing the CPB as a tool to 
generate single photons for state preparation and as a 
non-linear phase shifter to construct two qubit gates be- 
tween different molecular ensembles. 

For molecular ensembles to qualify as a quantum mem- 
ory we not only need fast read/ write operations but we 
also demand that the lifetime of an arbitrary ensemble 
superposition lipe) — a|Oe) + f3\le) is substantially longer 
than coherence times of solid state qubits. The lifetime 
of an ensemble qubits is limited by the single molecule 
decoherence time as well as mutual interactions. With 
expected single molecule decoherence rates of the order 
of a few Hz [l| we identify state dependent elastic and 
inelastic collisions [5] between molecules as the funda- 
mental limitation of the lifetime of a molecular ensemble 
quantum memory stored in a thermal gas of molecules. 

B. Dipolar crystal 

To avoid coUisional dephasing while still keeping the 
benefit from a collectively enhanced coupling we con- 
sider molecular ensemble qubits prepared in a high den- 
sity crystalline phase under ID or 2D trapping condi- 
tions. Dipolar crystals of molecules have recently been 
studied in Ref. ITj, where it has been shown that with 
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dipole moments aligned by an external DC electric field 
molecules are stabilized by repulsive dipole-dipole inter- 
actions Vdd — Cz/r^ in a plane perpendicular to the 
bias field. Attractive interactions along the remaining 
direction are suppressed by a strong transverse confine- 
ment [i3,[i8| . Under such conditions the low temperature 
physics of the molecules is characterized by the dimen- 
sionless parameter 



pot 



ao 



(4) 



which is the ratio between potential energy and kinetic 
energy for molecules of mass m for a given density 
n = l/(ao)'' and dimension d — 1,2. For 7^1 the 
dipolar repulsion wins over kinetic energy leading to the 
formation of a crystalline phase, i.e. small oscillations of 
the molecules around their equilibrium values. The for- 
mation of a dipolar crystal at high densities, i.e. in the 
limit where collisions are most damaging is in contrast to 
the familiar Wigner crystal for trapped ions or electrons 
where crystallization occurs at low densities. Numerical 
Monte Carlo simulations [13, [l^, [13, |21J have predicted 
a crystalline phase for 7 > 1 in ID and 7 > 20 in 2D. 
For typical experimental numbers a stable crystal of po- 
lar molecules is found for a lattice spacing oq of a few 
times 100 nm. 



C. Ensemble qubits in dipolar crystals 



III. ROTATIONAL ENSEMBLE QUBITS IN A 
HOMOGENEOUS DIPOLAR CRYSTAL 



In this section we consider the properties of ensemble 
qubits with N molecules prepared in a crystalline phase, 
and qubits encoded in collective rotational excitations. 
Our goal is to study the dynamics, and thus decoherence, 
of an initial rotational ensemble qubit \'ipe) — Q!|Oe)-f/?|le) 
under the influence of dipole-dipole interactions. We 
start with the simplest possible model of a homogeneous 
ID or 2D crystal, returning to questions of experimen- 
tal implementations and requirements (e.g. questions of 
transverse and longitudinal trapping potentials, and the 
choice of particular molecular states) at a later stage. 



A. Hamiltonian 

Let us consider an (infinite) homogeneous dipolar crys- 
tal of a given density corresponding to a lattice spac- 
ing flo, which is initially prepared in the qubit state 
|0e) — \gi ■ ■ -Qn)- We denote by r? the classical equilib- 
rium positions of the molecules, which form a linear chain 
in ID or a triangular lattice in 2D. As discussed in Sec. [Til 
the stability of the crystal requires 7 = Csm/h^ao 3> 1, 
where now C3 = ^'^/Atteq is determined by /Zg, the in- 
duced dipole moment of state \g). The dynamics of the 
system including internal and motional degrees of free- 
dom is given by the Hamiltonian 



Hmdc = J2(?^ + f^egh){e,\ ] + Vddi{r^}) , (5) 



In the following sections we address the question 
whether it is possible to achieve a stable molecular crys- 
tal and at the same time encode quantum information 
in ensemble superpositions |V'e) = Q;|Oe) + /3|le). In this 
context we distinguish between two types of ensemble 
qubits: (i) rotational qubits, i.e. the states |0e) and |le) 
introduced above, and (ii) spin qubits. In Sec. lIIII we first 
consider rotational ensemble qubits which are directly af- 
fected by state dependent dipole-dipole interactions and 
decay by phonon induced scattering processes out of the 
symmetric state |le). In Sec.[IV]we extend our model to 
molecules with an additional spin degree of freedom and 
study ensemble qubits encoded in collective excitations 
of two spin states \g) and |s) within the same rotational 
manifold. As spin degrees of freedom are essentially 
unaffected by dipole-dipole interactions, spin ensemble 
qubits in a MDC form indeed a highly protected quan- 
tum memory. However, a degrading of the spin ensemble 
quantum memory due to dipole-dipole interactions still 
occurs during gate operations when molecules are (virtu- 
ally) excited into the rotational state |e). To estimate the 
resulting gate fidelities under realistic experimental con- 
ditions we focus in Sec.|V]on a specific implementation of 
a ID dipolar crystal and include effects of an additional 
longitudinal confining potential into our model. 



with ID or 2D position and momentum operators de- 
noted by Vi and p,;, respectively, and the transition 
frequency between states |e) and \g). The Hamiltonian 
([5]) is the sum of the kinetic energies of the molecules, 
a bare molecular Hamiltonian for the internal (rota- 
tional) states, and the dipole-dipole interaction Vdd({ri}) 
which couples the internal and motional degrees of free- 
dom. With fi denoting the electric dipole operator of the 
molecule the dipole-dipole interaction is given by 



Vddi{v^}) 



87re, 



1 V- MrMj -3(n.,j -/i,)(nij -/x^-) 

2^ lr.--r.|3 ' 



where n^j = r^j / 



ra I is the unit vector pointing along 
r,; - r,-. 



direction — j., - ij 

To study the dynamics of ensemble states |0e) and |le) 
under the action of Hmbc we proceed as follows. Since in 
the crystalline phase molecules are located around equi- 
librium positions « r° we describe in a first step the ac- 
tion of the internal operator Vi = Vdd({ri }) on the qubit 
state |le). As the dipole-dipole interactions depend on 
the actual choice of rotational states \g) and |e) we start 
with a short overview on rotational spectroscopy of polar 
molecules in Sec. IIII 51 In Sec. IIII Cl we then discuss the 
action of dipole-dipole interactions on ensemble states 
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|0e) and |le). We find that at least in a homogeneous 
crystal the only effect of Vj is an energy shift for state 
|le) which does not destroy the decoherence of a qubit 
state \tpe)- Therefore, in Sec. IIIIEl we include molecular 
motion and write the position operator of each molecule 
as Ti = + Xi, with accounting for small fluctua- 
tions around the classical equilibrium positions r^. By 
expanding Vdd({ri}) i"^ powers of x^ we obtain the dom- 
inant contributions for the interactions between internal 
and external degrees of freedom which, for example, ac- 
count for state dependent forces on the molecules due 
to a difference of the induced dipole moments of states 
\g) and |e). The analysis of the resulting model given 
in Sec. IIII Fl finally allows us to estimate the ensemble 
quantum memory lifetime, Tg, for a wide range of sys- 
tem parameters. 

B. Rotational Spectroscopy 

Since the interactions between molecules depend on the 
actual choice of rotational states \g) and |e) we will first 
summarize the rotational spectroscopy of polar molecules 
in the presence of external electric fields 22]. To keep the 
discussion on a basic level we consider in this section only 
molecules like SrO or CsRb with a closed electron shell 
and a electronic ground state. In Sec. IIVI we extend 
our model to molecules with additional spin degrees of 
freedom. 

At sub-Kelvin temperatures with electronic and vibra- 
tional degrees of freedom frozen out the energy spectrum 
of a molecule is well described by the rigid rotor 
Hamiltonian Hm = -BN^ with N the angular momen- 
tum of the nuclei and B the rotational constant which 
is typically in the order of several GHz. In the presence 
of an external electric bias field E(,, polar molecules in- 
teract with the field via the dipole coupHng — /xEb, with 
fl the electric dipole operator of the molecule. In the 
following we choose our z-axis along the direction of the 
bias field, i.e. = Ebe^ and the total Hamiltonian is 
Hm = - fi.Eb. 

For a given electric field we label eigenstates of Hm 
as \N,M]y)Et with A/jv the eigenvalue of the operator 
Nz. In the field free case {Eb = 0) the eigenstates 
|iV, MAr)o are the usual angular momentum eigenfunc- 
tions and the corresponding anharmonic energy spec- 
trum. En = B{N + 1)N, with energy levels {2N + 1)- 
fold degenerate is shown in Fig. [2] a). In the pres- 
ence of an electric bias field {Ei, ^ 0) the dipole cou- 
pling mixes different angular momentum eigenfunctions 
and general eigenstates \N ^M^) E^, are superpositions of 
states I A^, Mn)o with different N but with the same Mn 
quantum number. The modified rotor spectrum as a 
function of the applied electric field Ei, is show in Fig. [2] 
b). The spectrum retains its anharmonicity over a wide 
range of electric fields values and only for E}, ^ B/^q 
with fiQ the axis fixed dipole moment of the molecule 
rotational excitations become approximately harmonic. 
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Figure 2: a) Rotational energy spectrum for zero bias field 
Eb = (solid lines) and the corresponding shifted energy 
levels at Ei, — SB/no (dashed lines), b) Energies spectrum 
and c) induced dipole moments /ii — {j|^z |i) for the lowest few 
eigenstates as a function of the applied bias field Et- Different 
curves for eigenstates \N, Mn) are labelled by (A'^, |Afjv|). 



The mixing of different rotational states |7V, Afjv)o in 
the presence of an electric field corresponds to an align- 
ment of the molecules along the field direction, and in 
general an eigenstate |0) = \N, Mn)^^ exhibits a fi- 
nite induced dipole moment {(j)\fj,z\4') along the z- 
direction. The magnitude and sign of the induced dipole 
moments depend on the specific state and the strength 
of the bias field Eb- The induced dipole moments for 
the lowest rotational states are plotted in Fig. [5] c) as a 
function of the electric field Ei,. We find that for an axis 
fixed dipole moment /io ~ 5 — 10 Debye and moderate 
electric fields, Ef, w 2B/^q, which typically corresponds 
to a few kV/cm, induced dipole moments of a few De- 
bye can be achieved. Note that in general induced dipole 
moments of non-degenerate states are different, although 
there are certain values of the bias field, so-called 'sweet 
spots' where induced dipole moments for two states 
are the same. 

The anharmonicity of the rotor spectrum for low and 
moderate electric fields allows us to choose two rotational 
states, e.g. \g) = \N,Mn)e, and |e) = \N',M'j^)e,, 
which are selectively coupled by a dipole transition to 
the fields of a microwave cavity of appropriate frequency 
and polarization. Selection rules for dipole transitions 
require AMn = 0,±1, while the additional restriction 
AA^ = ±1 applies only for vanishing Ei,. In the following 
the states \g) and |e) form our truncated single molecule 
basis which defines our ensemble states |0e) and |le), as 
introduced in Sec. HI 
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C. Dipole-dipole interactions 

So far we have discussed ensemble qubits on the level of 
independent molecules. As the crystalline phase is stabi- 
lized by mutual interactions between molecules, which for 
different internal states depend, e.g., on the different in- 
duced dipole moments, we proceed to study dipole-dipole 
interactions given in Eq. 

Consider the action of Vdd on the qubit states. For 
fixed positions we define the operator acting on the 
internal states, 

Aj(n,j) = fi. fij - 3(ny-/xJ(ny-/Xj.) . (7) 

We will simplify Dij{nij) under the assumption that its 
action is confined to the two molecule subspace Hij = 
{\ei,ej),ei = g,e} and that it is independent of 
when molecules are confined in the (x, ?;)-plane. This 
simplification is possible under the following conditions. 
First, due to the comparatively large separation between 
molecules of |ry| > 100 nm the dipole-dipole interaction 
is a weak perturbation on the scale of the rotational spec- 
trum (~ B). This is true for most choices of states \g) 
and |e) and allows us to restrict the action of the operator 
Dij{xiij) to the subspace Ti-ij- There are exceptions, e.g. 
the combination \g) = |iV, 0) and |e) = |7V',-|-1) where 
the degeneracy between the states |iV', 1) and |A^',— 1) 
would lead to resonant transitions out of Ti-ij- We can 
nevertheless include those combinations of states in our 
discussion when we assume that the degeneracy is lifted, 
e.g., by additional AC microwave fields. Our second as- 
sumption is that with molecular motion restricted to the 
(x, ?/)-plane, i.e., riij _L the operator Dij{nij) is in- 
dependent of riij. This condition is fulfilled for \g) and 
|e) being eigenstates of N^. For other states, e.g. |e) as 
eigenstate of the operator , the directional dependence 
of Dijip-ij) would lead to an additional x/y dependence 
in the models for 2D crystals derived below. For simplic- 
ity we ignore these cases in the following discussion. 

With these assumption we decompose Dij = Dlj + 

D^lj' ■ The first term, D^l^ , describes a state dependent 
interactions due to a difference in the induced dipole mo- 
ments, 6fi = — fj,g, with Hi = {i\nz\i). It is diagonal in 
the qubit basis and we can write it as 

D^f =i^^gh + S^l\e,){e,\){^lgl,+S^,\e,){e,\) . (8) 

^ (r) 

The second term, D^j' , accounts for resonant exchange 
of rotational quanta between molecules. Introducing a 
Pauli operators notation cr^ = \ei){gi\ etc., it is given by 

A?=D.(f7+a-+ar^+) , (9) 

with the matrix element D,, ~ 7]\{e\fjb\g)\'^ , where = 1 
for rotational states of equal quantum number Mj^ and 
7? = -1/2 for M'j^ = Mat ± 1. 



While the diagonal operator D\j acting on the qubit 
states leads to a state dependent energy shift, the oper- 

" (r) 

ator flips the rotational excitations between neigh- 
boring molecules. When motional degrees of freedom are 
included, both processes result in state dependent forces 
on the molecules. 



D. Effect of dipole-dipole interactions for fixed 
lattice positions 

Let us consider the simple situation where the 
molecules are frozen at lattice positions r". In this case 
the dipole-dipole interaction takes on the form 

In view of [V/,i?t] ^ both ensemble qubit states |0e) 
and |le) are eigenstates of Vj. Therefore, apart from a 
small energy shift, the internal part of the dipole-dipole 
interactions does not limit the lifetime of the ensemble 
qubit. This statement, of course, ignores inhomogeneity 
and finite size effects which depend on the specific ex- 
perimental setup (see Sec. |Vl. However, these imperfec- 
tions can in principle be avoided and do not constitute 
a fundamental restriction to ensemble quantum memo- 
ries in dipolar crystals. We conclude that in our model 
state dependent forces and the resulting entanglement 
with motional degrees of freedom is the primary source 
of decoherence. 



E. Effect of dipole-dipole interactions including 
motional couplings 

We return to the full Hamiltonian -ffiviDC given in 
Eq. (O which includes internal as well as external de- 
grees of freedom. 

1. Decomposition of the dipole-dipole interactions acting on 
internal and motional degrees of freedom 

With the assumption that the dipole-dipole interaction 
is confined to the subspace Hij = {\ei,ej),ei — g,e}, we 
write 

where we introduced the dimensionless operator Kij by 
= A^g (ly + K,j^ ■ (12) 

This decompositions separates Vdd into a part which is 
independent of the internal state and describes purely 
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repulsive interaction between molecules which stabilize 
the crystal. All state dependent properties of are 
contained in the operator Kij given by 

K^J = e (|e,) (e,| + |e,) (e, |) + k {a+aj + a+) . (13) 

Here e — {fie — fJ-g)/fJ-g the normalized difference of in- 
duced dipole moments, and k — Dr/fJ.^ the normalized 
coupling constant for resonant exchange processes. Note 
that in Eq. ()13|1 we have omitted the term |ei, e^) (e^, Cj \ 
which is negligible for a low number of rotational excita- 
tions, as is the case for our initial ensemble qubit. There- 
fore, for a given choice of states \g) and |e) we characterize 
dipole-dipole interactions by the induced dipole moment 
of the ground state fig, and the two dimcnsionless pa- 
rameters e and k. 

We rewrite the molecular position operators as = 
-I- Xi and expand Eq. pT|) in x^, so that the dipole- 
dipole interaction splits into three contributions, 

Vddi{r^}) = Vi + VEi{x^}) + Vi„t({x,}) , (14) 

with Vi (Ve) acting on internal (external) degrees of free- 
dom respectively, while Vint contains all remaining terms 
which couple the external and the internal dynamics. 

2. Excitons: rotational excitations hopping on the lattice 

Let us first return to the internal operator Vi defined 
in Eq. (jlOp , where by neglecting a global energy shift we 
can replace operators Dij by fx^Kij. The operator Kij 
given in Eq. (|13p preserves the total number of molecules 
in state |e), but allows a propagation of rotational exci- 
tations on the lattice. In the limit of a low number of 
rotational excitations we diagonalize Vj by introducing a 
set of collective operators defined by 

4|0e) ^ |lk) = 1/ViV^ e"^"-?|5i ... e, ... 5^), 
j 

with a wave vector k restricted to the first Brillouin zone 
of the lattice. The operators i?k fulfill (approximate) 
bosonic commutation relations, [RiujRI^,] — ^kk'j and in 
the following we refer to states created by i?]^. as 'ex- 
citons'. This nomenclature is based on the similarities 
of rotational excitations with localized Frenkel excitions 
in organic crystals [1^. We then identify our qubit 
state |le) as associated with the zero momentum exci- 
ton, RI = RI^q- Including the energy offset fwjeg the 
dynamics of these excitons is given by the Hamiltonian 

i?exc = J2 f^eg\e^){e,\ + Vi ^ ^(k)4^k ■ (15) 

i k 

The energy band of rotational excitations, i?(k), is given 
in Eq. ([2T|) of Sec. IIIIFl where we will discuss it in more 
detail. For the moment we simply note that Hamiltonian 
Hexc is diagonal in k, such that [H cxc, Rl] — E{0)Rl as 
already pointed out at end of Sec. IIII CI 



3. Phonons 

In a next step we consider the dynamics of exter- 
nal degrees of freedom of the molecules which is deter- 
mined by the interaction VE{{x.i\) and the kinetic en- 
ergy i?kin = ^i'P'i I'^'iTL. Since the first order expansion 
of VB({xj}) vanishes due to the definition of equilibrium 
positions r°, the first non- vanishing contribution is of 
second order in x^ and given by 

V (f^ l\ \" 5[(x,;^Xj-)TlO.]'-[x,-X,f 

(16) 

The quadratic interaction between molecules causes col- 
lective oscillations (phonons) in the crystal described by 

the Hamiltonian iJphon = ^^kin + VE{{^i})- As i/phon 
is simply a set of coupled harmonic oscillators it can be 
written in diagonal form 

i?phon = ^A(q)4(q)aA(q) . (17) 
q,A 

Here wc introduced the annihilation (creation) operators 
a\{q) (a|(q)) for phonons of quasi momentum q and fre- 
quency u!\{q). In 2D, the index A labels the two different 
phonon branches. The phonon modes in the dipolar crys- 
tal are acoustic phonons. A discussion of the frequency 
spectrum is given in Sec. IIIIFl 



4- Exciton-Phonon Interactions 

The remaining terms of Vdd({ri}) which can not be 
decomposed into purely internal or external operators are 
summarized in 14nt({xi}). The first non- vanishing order 
of 14nt({xi}) is linear in the operators x^ and is given by 

^"'(^'''^^^ ^ " 8^ ^ WZA^ (^'^ - ® ■ (18) 

It describes a state dependent force on the molecules and 
entangles internal and external degrees of freedom. In the 
following we introduce a new symbol iJjnt = Vint({xi}) 
and rewrite Eq. (|18|) in terms of exciton operators i?k 
and the phonon operators Oq. We obtain an interaction 
Hamiltonian of the form 

Hn,t = ^^^(q' k)[«A(q) + aU-q)] , (19) 

k,q,A 

which describes scattering processes from state |k) into 
state |k-|-q) under the absorbtion (emission) of a phonon 
of quasi momentum q (— q). We postpone a discus- 
sion of the explicit form of the coupling matrix elements 
MA(q,k) to Sec. HUB 
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5. Summary 

In summary, we have shown that the dynamics of a 
molecular dipolar crystal given by Hmuc iii Eq. ([S]) con- 
tains the three contributions, 



with 



-ffMDC — Hcxc + Hphon + -^int , 
k 

Hphon = ^hajxiq)a{{q)ax{(i), (20) 
Hint = ^ AfA(q,k)[aA(q)+a^(-q)]4+qi?k, 

k,q,A 

which is the (minimal) model which describes the evo- 
lution of ensemble qubits in a self-assembled molecular 
dipolar crystal. While explicit expressions for the energy 
dispersion -B(k), the phonon spectrum ijJ>(q) and the 
scattering matrix elements M\{q, k) are give in Sec. lIIIFl 
for the ID and 2D crystal we first note the general struc- 
ture of -ffiviDC- The ensemble operator i?| is an eigcnop- 
erator of i?oxc as well as i?phon, and therefore, apart form 
an energy shift, the first two lines of Eq. (|20|) preserve the 
coherence of a qubit superposition [ipe) = a\Oe) + /3|le)- 
The third line of Eq. (|20|) , iJint , leads to phonon assisted 
transitions from the symmetric qubit state |le) = |k = 0) 
into orthogonal states |k ^ 0). This loss process is the 
dominant source of decoherence for a qubit state I^Ae) and 
in Sec. lIII Gl we calculate the resulting lifetime for the 
ensemble quantum memory. 



F. Molecular dipolar crystals in ID &: 2D 

In this section we discuss the exciton dispersion £'(k), 
the phonon spectrum lox (q) and the coupling matrix ele- 
ments M\{q, k) which determine the properties of Hamil- 
tonian i?MDC given in Eq. (j20|) . For a lattice spacing oq 
we express those quantities in terms of the dipole-dipole 
energy Udd = A*g/(47''eoag) and the dimensionless param- 
eters 7, e and k. Our focus is placed on the ID crys- 
tal where we derive analytic expressions for the relevant 
quantities. For the 2D crystal we present numerical re- 
sults and identify the main differences compared to the 
ID case. The derivations of the following results can be 
found in App. [Xland the main results are summarized in 
Fig. [3] (ID) and Fig. g] (2D). 



1. Excitons 

The energy spectrum of excitons in the dipolar crystal, 
i?(k), contains three contributions. 



where the dimensionless band structure J(k) is defined 
in App. |X] in Eq. (|A3p . While a finite difference in the 
dipolc moments, e 0, only causes a shift of the tran- 
sition frequency with J(0) — 2^(3) in the ID case and 
J(0) ~ 11.034 in 2D, the resonant exchange processes 
proportional to k lead to the formation of an exciton 
band structure as shown in Figs. [3] and 01 For the ID 
crystal the explicit expression for J(fc) is given in App. |^ 
in Eq. (IA4p and we find that in the long wavelength limit 
fc ^ it exhibits a non-analytic behavior. 



J(fc) ~ J(0) -3/2[l - 3/21og(fcao)](fcao)^ 



(22) 



which is a consequence of the slow decay of dipole-dipole 
interactions. For the 2D crystal the long-range char- 
acter of dipole-dipole interactions is even more appar- 
ent and results in a linear dispersion ECk) — E{0) ~ |k| 
for small |k|. The total width of the energy band is 
AE = 7C(3)/2 X \K\Udd in ID and AE ~ 13.37 x \K\Udd 
in 2D. Note that for rotational states with k positive the 
band structure is 'inverted' and the long-wavelength ex- 
citations have the highest energy. 



2. Phonon spectrum 

As shown in App. |^ the spectrum of the acoustic 
phonon modes in the self-assembled dipolar crystal is of 
the general from 



fkj xiq) = \/-Uddfx{q) . 



(23) 



In the ID crystal the dimensionless function f{q) defined 
in Eq. (|A8[) has a long wavelength limit f{q) ~ •\/l2C(3)x 
(jflo, and a maximum value of /(tt) = V93C(5)/2 ~ 6.94. 
The phonon spectrum is therefore characterized by the 
sound velocity c = ^12^(3)/7 aoUdd/fi and the Debye 
frequency, ljd = ^(^)i 



hujn — 



/93C(5) 



27 



X Udd- 



(24) 



The full phonon spectrum is plotted in Fig.[3]b). In the 
2D crystal there are two acoustic phonon branches, A = 
1,2 and the corresponding dimensionless spectra /A(q) 
are plotted in Fig. [4]b). In 2D the maximum phonon 
frequency is Huid — 8.22 x Udd/y/^- 



3. Exciton-phonon interactions 

Excitons and phonons interact via iJint given in the 
second line of Eq. ([20)) . We write the coupling matrix 
element as 



MA(q,k): 



1 



E{k) = hweg + Udd[eJ{Q) + «;J(k)] , 



(21) 



.Udd 



iegx{q) + K[9x{k+q)~gxim ^ 
(25) 
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Figure 3; Dipolar crystal in ID: a) Dimensionless band struc- 
ture J{k) (see Eq. (I21|l '). b) Dimensionless phonon spec- 
trum /(g) (see Eq. 123}). c) Function g{q) which enters in 
the expression of coupling matrix elements M [q, k) given in 
Eq. 




|k|ao M K 



where we introduced an additional dimensionless func- 
tion gx{cO defined in App. E]in Eq. \kl2\ . For the ID 
crystal the explicit expression of g{q) is given in Eq. (|A13|) 
and plotted in Fig. [3] c) . The matrix element M\ (q. k) 
contains two contributions. The first is proportional e 
and describes a phonon induced (on-site) energy shift of 
a molecule in state |e) due to a difference in the induced 
dipole moments. This type of interaction does only de- 
pend on the transferred momentum q and is familiar from 
polaron models discussed in solid state physics [l^l . The 
second contribution proportional to k describes phonon 
induced hopping of excitons. The coupling matrix ele- 
ments for this process also depend on the initial exciton 
state |k). In the long wavelength limit, |q| 0, both con- 
tributions scales as ~ \/|qJ, and scattering events with 
low momentum transfer are suppressed. 

Without going into the details of AIx{q^, k) wc point out 
two properties which are relevant for the discussion be- 
low. First, for long wavelength excitons the total strength 
of the exciton phonon interaction is in the order of 

|k| ^ : O(iJint) - k + n\UM ' , (26) 

and therefore only weakly suppressed by the parameter 
7. This means that in general the exciton-phonon in- 
teraction has a considerable effect on the dynamics of 
a molecular dipolar crystal even deep in the crystalline 
phase with 7^1. However, for a specific choice of ro- 
tational states and values of the electric bias field where 
e-|-K = is fulfilled, long wave length excitons completely 
decouple from the phonon modes. We come back to this 
point at the end of Sec. IIII Gl 

G. Lifetime of rotational ensemble qubits 

Based on the structure of Hmuc given in Eq. (PO)) and 
the discussion of its properties in subsection IIII Fl we 



Figure 4: Dipolar crystal in 2D: a) Dimensionless band struc- 
ture J(k) (see Eq. (HI}), b) Dimensionless phonon spectrum 
/A(q) (see Eq. ((23}). Band structure and phonon spectrum 
are plotted for k vectors along the symmetry lines of the first 
Brilluoin zone with symmetry points indicated in the inset of 
plot a). 

finally return to the original question of the lifetime of an 
ensemble qubit prepared in a state |^e) = ct\Oe) + /3|le)- 
As already mentioned above we find that the dominant 
decoherence mechanism arises form the decay of the state 
|le) due to phonon assisted scattering into orthogonal 
states |k 7^ 0). This means we can identify the lifetime 
of the ensemble quantum memory, T^, with the lifetime 
of the state |le). 



1. Ensemble qubit decay 

In the following we consider the situation where at 
some initial time t = the system is prepared in the 
state po = |0e)(0e| (8) pTwith px the equilibrium density 
operator of the phonon modes for a temperature T. At 
time i = 0+ we instantaneously excite the molecules into 
the qubit state |le) and calculate probability Pe{t) to find 
the system in state |le) after a waiting time t. As we are 
only interested in times t where Pe{t) is still close to 1, 
we can use second order perturbation theory and obtain 

P,(t)~l-A J' dt' J' dTY,\Mx{q,0)\^X 

[(iV(u;A(q)) + l)cos(!l-(q)T) +iV(wA(q))cos(fi+(q)T)] . 

(27) 

Here N{uj) ~ l/[exp{hLu/kBT) — 1] is the thermal occu- 
pation number for phonons of frequency ui and fl^ (q) = 
[£■(0) - E{(i)]/h ± tj(q). For very short times Eq. ^ 
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leads to a quadratic decay of the excited state probability, 
Pe(t) =i 1 - T^^^^ (28) 
with a characteristic rate W defined by 

= l^'^^(q'0)l'(2^(^^(q)) + 1) ■ (29) 

For long times the decay of Pe (t) turns into a linear func- 
tion of t, 

FeW^l-rt, (30) 
with the decay rate F given by Fermi's Golden Rule, 

r = § E 0)1' [(^K(q)) + l)Smq) c.,(q)) 

(31) 

The crossover time tc between the quadratic and the lin- 
ear regime is roughly given by (ic)~^ ~ max.{AE/h, o^^i}, 
with the exciton band width AE and the phonon Debye 
frequency ujd discussed in Sec. IIII Fl As long as the as- 
sociated decay probability Pc — W'^tj. is much smaller 
than one, the application of Fermi's Golden Rule is valid 
and we obtain a lifetime Te — l/F. We refer to this case 
as the weak coupling regime. Otherwise, for Pc ~ 1 or 
in the strong coupling regime, the qubit decay is deter- 
mined by the quadratic formula given in Eq. (|28p and 
we identify Te = 1/W. At low temperatures we find 
that when the exciton band width AE is larger than 
ftWD, Pc ~ (e + i^Y I {i^'^ while in the opposite case 
Pe^{e + nf. 



2. Quadratic decay: strong coupling regime 

In the strong coupling regime the initial non-energy 
conserving transitions out of the qubit state |le) already 
lead to a strong quadratic reduction of Pe{t)- The cor- 
responding rate W defined in Eq. ([M]) can be written 
as 



W 



2, _ulA^ + >^? 



Udd 



(32) 



with d = 1,2 the dimension of the crystal and Td{T) a 
numerical integral 



d^q |5A(q)|' 



bzVbz fx{q) \ef^(^')/^-l 



1 



(33) 

In the two limiting cases this integral behaves as Id(T — )■ 
oo) = 0(1) andX(;(T ^ 0) = 0{t). Neglecting numerical 
constants arising from the exact evaluation of TdiT) we 
can summarize the estimated qubit lifetime as 

(Te)-i « E^\l±ji X max{l,7VWt^<id}. (34) 



3. Fermi's Golden Rule: weak coupling regime 

In the weak coupling regime only energy conserved 
transitions with a certain transferred momentum qo lead 
to the decay of Pe{t). According to Eq. (pij) qo is defined 
by r^(qo)±a;A(qo) = 0, where the positive (negative) sign 
applies for a, k < (k>0) corresponding to phonon ab- 
sorption (emission) processes. Based on the discussion 
on the exciton and phonon spectra given in Sec. IIII Fl we 
distinguish between two cases. 

In the first case the resonant k-vector is different from 
zero, qo 7^ 0. Such a situation occurs ID for k/^/j^ 1, 
when the exciton bandwidth exceeds the phonon band 
width, AE > fvjjD- In 2D, due to the linear long- 
wavelength limit of E{k) this situation occurs only for 
a negligible parameter regime where AE ~ fiWD. The 
decay rate for ID is 



Udd (e + ^)^ 

h 



X C{qo) X [N{oj{qo)) + Q{k)] , (35) 



with C(go) = 2g^iqo)/(.f{qo)\J'{qo) + fiqo)/Vl^\) ^ «• 
The Heaviside function Q{k) takes into account that only 
for K > there is a finite decay rate at T = 0. 

In the second case the resonance condition is only ful- 
filled for qo = 0. This situation is in general true for 
a 2D crystal and in ID for k/^/j < 1. The decay rate 
then depends on the |q| — > limit of the summand in 
Eq. dSIl)- As |M(q, 0)p ~ |q| vanishes in the long wave- 
length limit, we only obtain a finite value for F in ID and 
for a non-zero temperature. The resulting rate is 



ID 



(e + K) 



(36) 



In 2D, due to the additional factor of |q| in the density 
of states, we obtain T2D = 0. 



4- Discussion: rotational ensemble quhits 

In summary we have found that the lifetime Te of a ro- 
tational ensemble quantum memory in a molecular dipo- 
lar crystal is primarily determined by the dipole-dipole 
energy Udd = Mg/(47i'eoao) and depends further on the 
relation between the dimensionless parameters 7, e and 

K. 

For low temperatures the dipole-dipole energy Udd is 
limited from below by the condition Udd > 7cfi.^/C|m'^, 
with 7c the minimum value of 7 which guarantees the 
stability of the crystal. For an upper bound we use the 
critical value 7c = 20 predicted in Ref. [l3, [Mj for the 2D 
crystal. For the nominal parameters fig = I Debye and 
m — 100 amu we obtain Udd/h > 260 kHz. As this bound 
scales like ^ 1/fig choosing molecular states with larger 
induced dipole moments fig ^ 5 D the scale of F can in 
principle be reduced below the kHz regime. However, for 
a system at finite temperature Udd is also bound from 
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\N, Mn) 




Ma bo] 


K 


e 


a) 


\g) = 11,0), |e) = |2,0) 


3.05 


-0.16 


10.5 





b) 


\g) = |l,0),le) = 13,0) 


3.91 


0.09 


1 





c) 


\g) = 10,0),le) = 11,0) 


8 


0.75 


0.1 


-0.7 


d) 


Ig) = 10,0),|e) = |1,1) 


5 


0.68 


-0.24 


-0.29 


e) 


|(7) = |0,0),|e) = |l,0) 


1.44 


0.39 


1.51 


-1.51 


f) 


\g) = |l,0),|e) = |3,0) 


3.44 


-0.13 


0.39 


-0.39 



Table I: Different combination of rotational states \g) and 
|e) and the corresponding values of e = (/ie — M9)/M9 ^iid 
K = Dr/y?g. States are specified by the quantum numbers 
|A'', Mat) and the value of the external bias field, Et,. 



below by fcsT. For a 2D crystal the melting temperature 
of a dipolar crystal is Tm ~ 0.09Udd/kB [25|, while in ID 
we expect a stable crystal for ksT < Udd (see Sec. IV)) . 

For a given Udd the lifetime of the ensemble qubit de- 
pends further on the dimensionless parameters k and e 
which in turn depend on the rotational basis states \g) 
and |e) and the value of the applied bias field Ej,. In 
Table U we have listed a few specific choices of rotational 
states \g) and |e) and the corresponding values of k and 
e. For certain states we find a 'sweet spot' of the elec- 
tric bias field where the induced dipole moments of two 
rotational states are the same and e vanishes. The value 
of K can be significantly reduced by choosing states \g) 
and |e) with a (to a high degree) dipole forbidden tran- 
sition, e.g. N' = N ± 2, with the cost of also lowering 
the coupling to the microwave cavity. An ideal situation 
occurs when both e and k are large but exactly cancel 
each other, e.g. e -I- k = 0. For such a 'magic' configu- 
ration the dipole-dipole interactions of the triplet states 
\gg) and \ge) + \eg) are exactly the same, which makes a 
symmetric excitation jig) insensitive to phonon-induced 
fluctuations of dipole-dipole interactions. Examples for 
such phonon-decoupled states are given in row e) and f) 
of Table [D An additional interesting observation is the 
absence of energy conserving transitions in a 2D crystal, 
which in the weak coupling regime implies T2D = 0. A 
detailed investigation of this case and corrections due to 
higher order terms in the interactions Hamiltonian 
are the subject of future research. 

In conclusion we find that the phonon induced decay 
of rotational ensemble qubits can cover a wide range of 
values and is in the end to a large extend determined by 
the specific experimental implementation of this system. 
When trapping of molecules can be achieved independent 
of the rotational state, e.g. in optical or magnetic traps, 
we can find certain magic configurations in which rota- 
tional ensemble qubits are completely decoupled from the 
phonon bath. In other cases, e.g. for electrostatic traps, 
we are limited to a certain set of trappable rotational 
states and typically we are not able to fulfill the decou- 
pling condition e -|- k = 0. In such a situation it is not a 
good choice to encode quantum information in the rota- 
tional degree of freedom. 



IV. SPIN ENSEMBLE QUBITS: 

In this section we extend our discussion of molecu- 
lar ensemble qubits to molecules with spin or hyperfine 
states. This additional internal degree of freedom which 
is not directly affected by electric dipole-dipole interac- 
tions allows us to overcome some of the limitations of 
purely rotational qubits. In particular, it is then possible 
to use the rotational degree of freedom to trap molecules 
in electrostatic potentials and encode quantum informa- 
tion in different spin states, or in contrast to use spin 
states for magnetic trapping and encode quantum infor- 
mation in rotational degrees of freedom. 

A. Rotational spectroscopy of molecules 

In general the rotational spectroscopy of polar 
molecules involves spin-rotation interactions between the 
rotation of the nuclei N and the spin of unpaired elec- 
trons S as well as hyperfine interactions between S and 
the nuclear spin I |22l |. However, to keep the following 
discussion on a basic level we consider in this paper only 
the case of a molecules with a single unpaired electron 
with spin S — 1/2 and no nuclear spin / = for which 
hyperfine interactions are absent. As, to our knowledge, 
non of the polar molecules studied in current experiments 
falls into this category [2^ we explain below how the ar- 
guments presented for the molecule can be applied for 
more complicated molecules with / > as long as the hy- 
perfine interactions are small compared to the rotational 
constant B. 

In the presence of a bias field = Ei,ez the Hamil- 
tonian for a molecule in the vibrational ground state 
is 

HM = BN^+^srNS-fi,Eb, (37) 

with the spin- rotation coupling constant, typically in 
the order of 100 MHz. For a vanishing bias field the spin- 
rotation coupling lifts the degeneracy of | A^, Mn)o states 
and new eigenstates of Hm are given by |A^, S; J, Mj)o 
with J = N -f S the total angular momentum. For mod- 
erate and strong electric fields Ei, > B/hq the dipole 
coupling typically exceeds the spin-rotation coupling and 
eigenstates of Hm are approximately given by product 
states \N,MN)Et \^s)- Therefore, to a good approxi- 
mation energy splitting and induced dipole moments are 
determined by the rotational component only and we re- 
fer to the discussion given in Sec. IIII Bl Especially, for 
states with Mn = different spin components are degen- 
erate and have exactly the same induced dipole moments. 
For other values of Mn 7^ different spin components are 
split by the spin rotation coupling, jsr, but still have to 
a good approximation the same induced dipole moments. 
As an example we plot in Fig. [5] the spectrum of Hm at 
the 'sweet spot' of the bias field, E^ = Es = 3.055/^0- 

For molecules with a finite nuclear spin / > the 
physical picture is quite similar although the resulting 
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Figure 5: Spectrum of a molecule for a bias field Eb = 
3.05-B//j,o. Eigenstates which are approximately of the form 
|iV, Mjv)i3t ® \Ms) are individually labelled by {Mn,Ms) 
while the quantum number TV is indicated only for the corre- 
sponding manifold. 



spectrum is more involved. In the limit of 7.,^ — > the 
hyperfine interaction Hamihonian H^f couples the elec- 
tron spin and the nuclear spin to a combined angular 
momentum F3 = S + I and eigenstates are of the form 
\N,Mm)e,, ® \FztMp.^). For a strong electric bias field 
and for Mjv = this factorized form is approximately 
conserved even for finite ^sr as the spin-rotation coupling 
is too weak to mix different rotational states. For non- 
zero the diagonal part of the spin-rotation coupling, 
^srNzSz, does in addition lead to a mixing of different 
states iF^jMps)- This general picture of the spectrum 
holds as long as \Hhf\,jsr B < noEt while details 
depend on the exact relation between hyperfine and spin 
rotation coupling. 



B. Spin ensemble qubits in MDC: a protected 
quantum memory 

The additional spin degree of freedom allows us to en- 
code quantum information in ensemble qubits |0e) and 
|le), where the corresponding molecular basis states \g) 
and |e) have different spin components. In the following 
we explain in two examples how the additional spin de- 
gree of freedom can compared to purely rotational states 
improve the lifetime of molecular ensemble qubits in a 
dipolar crystal. 



1. Long-lived ensemble quantum memory 

We first consider molecular basis states \g) — Ei,® 
\Ms = -1/2) and |e) = |2, 0) ® Ws = 1/2) for the bias 
field Eh = Es- As we discuss in Sec. |V] this is an ex- 



ample for two states which can be both trapped in elec- 
trostatic fields. The states \g) and |e) have different ro- 
tational components and different spin components. At 
this point it is important to note that the factorization 
of eigenstates of Hm in a rotational and a spin com- 
ponent is only approximate. The spin-rotation coupling 
still admixes a small fraction of the other spin component 
which results in a finite transition dipole matrix element 
I (el/xjg) I, which actually allows us to couple states \g) and 
|e) to the microwave cavity. The dimensionless transition 
matrix element 0^: = |(e|/ia;|(7)|//io which determines the 
coupling to the cavity as well as the parameter /t depends 
on the ratio ^sr/B. For the typical case of ^sr ^ B 
we find in numerical calculations Qx — 2.5 x ^g^/B. 
For molecules like CaBr with ^sr/B — 0.03 the spin 
forbidden transition is roughly a factor 7 smaller than 
the corresponding spin-allowed transition. In contrast 
the dimensionless parameter k ~ 0^ is reduced by 
a factor {'jsr/B)'^ and exciton phonon interactions are 
strongly suppressed. For our specific example we obtain 
K ~ 250 X {'jsr/B)^ <C 1, which should be compared to 
K « 10 for spin conserving transitions (see row a) in Ta- 
ble |T|. Therefore, the spin degree of freedom provides an 
additional knob to change the interaction parameters k 
and e independently and tune the crystal into the weak 
coupling regime where decay rates are highly reduced (see 
discussion given in Sec. IIII G]) . 



2. Protected spin ensemble quantum memory 

To improve the lifetime of the ensemble quantum mem- 
ory even further we consider in a next step the states 
\g) = \1,0)eMMs = ~1/2) and \s) - |1, O)^,® |Ms = 1/2). 
As l^) and \s) have the same rotational component we 
have introduced the new notation \s) to distinguish this 
state from the rotationally excited state |e) discussed 
above. While the states \g) and \s) can not directly 
be coupled with a single microwave photon, we can still 
achieve a coupling to the microwave cavity by a two pho- 
ton process involving an additional classical microwave 
field. 

As explained above the two spin states \g) and |s) 
have the same induced dipole moment, e = 0, and in 
addition there are no resonant dipolar exchange pro- 
cesses between molecules in states \g) and \s) which im- 
plies that also k — 0. Therefore, quantum informa- 
tion encoded in spin ensemble qubits |0s) = |0e) and 

|ls) = Sl\Os) = l/VNj2i\9i---Si---9N) is naturally 
protected from dipole-dipole interactions and the result- 
ing phonon induced decay. Higher order spin flip pro- 
cesses due to virtual excitations into higher rotational 
states H^, il^ , which are not included in our model scale 
as ^ ^Q/(167r^eQao) x ^l^/B^ and even for ao = 50 nm, 
this rate is only in the order of a few Hz. Therefore, the 
lifetime of spin ensemble qubits is mainly limited by the 
dephasing rate of the spin (or hyperfine states) of a single 
molecule. Similar to the case of cold atoms [29, wOl or 
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trapped ions [3l|, [s^l dephasing rates of molecular hyper- 
fine states of below 1 Hz should be achievable. 

While highly protected from dipolar interactions spin 
ensemble qubits are not directly coupled to microwave 
photons and the storage and retrieve operations require 
a two step process 

(a|0,) +/3|lc)) ^ {a\Oe) + P\le)) ^ (a|0,) + /3|1,)) , 

which involves the rotationally excited ensemble state 
|le). During this process the rotational ensemble qubit is 
affected by interactions with phonons which cause a de- 
cay into orthogonal states |lk) as discussed in Sec. IIII G\ 
Therefore, the overall fidelity of an ensemble quantum 
memory in a MDC is still affected by exciton-phonon in- 
teractions. However, by employing a Raman process as 
discussed in Ref. rotational states are only virtually 
populated and interactions with phonons are suppressed 
by the detuning A — lUc ^ LOeg- As we have not yet ana- 
lyzed the details of the cavity-ensemble coupling we post- 
pone the discussion of gate fidehties to Sec. IV Dl where we 
study swap operations between microwave photons and 
spin ensemble qubits for a specific setup. 

V. DIPOLAR CRYSTALS IN A TRAP: 
INTERFACING MOLECULAR ENSEMBLES AND 
CIRCUIT QED 

In Sec. mil and Sec. IIVI we have studied an idealized 
homogeneous MDC and identified the exciton-phonon in- 
teraction as the main source of decoherence for ensemble 
qubits. We have shown that for certain choices of molec- 
ular states, in particular states with a different spin com- 
ponent, phonon induced decay processes are suppressed 
and a highly protected ensemble quantum memory can 
be realized with this system. However, so far we have 
ignored questions related to the experimental implemen- 
tation of a MDC, especially questions of transverse and 
longitudinal trapping potentials. In this section we study 
the properties of molecular dipolar crystals under realis- 
tic experimental conditions, especially in the presence of 
an additional longitudinal trapping potential. 

For a coherent integration of a MDCs with a cir- 
cuit QED system as shown in Fig. [T] trapping of polar 
molecules close above the chip surface must not affect 
coherence properties of the superconducting device. In 
particular trapping techniques which require the applica- 
tion of strong magnetic fields or intense laser fields raises 
compatibility questions with high Q-values of supercon- 
ducting strip line cavities [lB|. Therefore, to achieve a 
strong transverse confinement which is compatible with 
long coherence times of the microwave cavity, we focus in 
this paper on a scenario where molecules are trapped by 
an electrostatic potential as recently proposed in Ref. [l| . 
In Sec. IV Al with briefiy outline the basic idea of elec- 
trostatic trapping of polar molecules and show that this 
specific trap design will restrict our choice of molecu- 
lar basis states to a very limited set of rotational states. 



In Sec. IV BI we then study the properties of excitons and 
phonons in a quasi ID trapping configuration with an ad- 
ditional harmonic confinement potential along the crystal 
axis. As the confining potential removes the translational 
symmetry of the crystal, symmetric excitations |le) are 
no longer eigenstates of i/cxc which opens a new decay 
channel for rotational ensemble qubits. In Sec. IV CI we 
use the spectra of longitudinal and transverse phonon 
modes to discuss the stability of a ID crystal in the 
case of finite temperature and a finite transverse trap- 
ping frequency. Finally, in Sec. IVE)] we use these results 
to discuss state transfer fidelities between a microwave 
cavity and spin ensemble quibts under realistic experi- 
mental conditions. 



A. Electrostatic confinement of polar molecules 

In the following we consider a molecule as discussed 
in Sec. IIV Al in the presence of a bias field Efc(r) which 
now depends on the position of the molecule. The Hamil- 
tonian is 

HM{r) = + 7,,NS - /xEb(r) . (38) 

To achieve trapping wc consider an electric field of the 
form Efc(r) — {Eotf + Et{r))ez, with a large offset field 
Eos and a small trapping field Et{r) with < Et{r) <C 
Eos and Et{r) = at the center of the trap. Without 
going into the details of the actual trap design we here 
envision an elongated version of the electric z-trap trap 
proposed in Ref. [1| which would produce an electric field 
configuration of approximately this from. 

As long as the position dependent trapping field 
Et{r) is small compared to the offset field Eos we can 
use a Born-Oppenheimer argument and write Hamilto- 
nian ([55]) as 

HM{r) ~ + Vt,n{r))\n){n\ , (39) 

n 

with |n) (En) denoting eigenstates (eigenvalues) of 
Hm for the bias field E;, = EosBz and V(,„(r) — 
— {n\iJ,z\'n)Et{r). As the electric field has a local min- 
imum at the trap center only 'weak field seekers' with 
Mn ^ < are trapped in this potential. From 

the discussion in Sec. IIII BI we find that this restriction 
limits our choice of molecular basis states to states with 
Mn = and moderate electric fields. 

As indicated by the index n the trapping potential Vt.„ 
depends in general on the internal eigenstate \n). In the 
following we avoid this dependence by tuning the offset 
field to the 'sweet spot', Eos = Es- We choose the two 
spin states in the N — 1 manifold, jg) = 1 1, 0)^;^ (8) \Ms — 
-1/2) and \s) = \1,0)es <^ l^s = +1/2) as our single 
molecule basis states for spin ensemble qubits. Employ- 
ing a two photon process the two spin states can be cou- 
pled to the microwave cavity via a third, rotationally 
excited state |e) = |2,0)es ® l^s = -1/2)- Restricted 
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to those three basis states we can write the molecular 
Hamiltonian as 



HMiv)^hLO,g\e){e\ + Vt{r)^ 



(40) 



with Vt (r) a state independent trapping potential for the 
molecule. Below we consider a quasi ID trapping config- 
uration with 



1 



Vt{r) = —rnvx' + v 



(41) 



Here v^^ is the trapping frequency for the strong transver- 
sal confinement and v the trapping frequency for 
an additional weak confinement along the crystal axis. In 
this electrostatic trap the transverse trapping frequencies 
can be as high as :^_L/27r « 1 — 10 MHz [l|]. 



B. MDC under quasi ID trapping conditions 

In this section we study the properties of a MDC in a 
quasi ID trapping configuration where compared to the 
discussion given for a homogeneous system in Sec. Illll we 
add a finite longitudinal and transverse trapping poten- 
tial Vt(r) as given in Eq. (|^T|l . Restricted to the basis 
states |f;), \s) and |e) identified in Sec. IV Al we extend 
the crystal Hamiltonian i?MDC given in Eq. ([5]) by the 
trapping potential Vt (r) and write the total Hamiltonian 
for the inhomogeneous MDC as 



MDC 



E 



2m 



+ Vt{Y.j) + nWeg\e^){e^\j + Vdd{{Yi\) , 

(42) 

with the dipole-dipole interaction Vdd{{'Yi\) given in 
Eq. ([6]). We proceed as in Sec. IIIII and describe dipole- 
dipole interaction by the induced dipole moment /Xg and 
the dimensionless operator Kij. At the 'sweet spot' with 
e = the latter is given by 



K,. 



^^{\9^.e:j){e^,g■j\ + \e^,g;j){g^,eJ\) . (43) 



For our specific choice of rotational states \g) and |e) we 
find K ~ 10.5 (see Table U example a)). Note that we 
can omit resonant dipolar interaction between states \s) 
and |e) as long as most molecules remain in state \g). 

Assuming a crystalline phase we replace the molecular 
position operators by ~ r° + x^, with r° — (a;°,0,0) 
and a;^ the classical equilibrium position along the crys- 
tal axis. As molecules are confined by an additional 
longitudinal trapping potential the positions are no 
longer equidistant. In contrast to the discussion given 
in Sec. IIIII we here also include fiuctuations along the 
transverse directions, e.g. = (xi,yi, Zi), to study the 
effect of a strong but finite transverse confinement. Ex- 
panding Hmuc up to the lowest relevant order in x; and 
decompose the crystal Hamiltonian into an exciton part, 
a phonon part and exciton-phonon interactions, 



MDC 



phoii 



(44) 



In terms of equilibrium positions the exciton Hamil- 
tonian is given by 



huj 



\ei){ei 



87reo ^ 



,.013 



(45) 



The phonon Hamiltonian contains now both longitudinal 

as well as transversal phonons, ffphon = ^^phon + ^jAon- 
Primarily we are interested in the longitudinal part which 
is of the form 



phon 



2 

2m 

3 
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87reo 



E 



{Xj Xj ) 
|„0 _ ^015 



(46) 



Although our focus is placed on the quasi one dimensional 
regime with transverse motion frozen out we include for 
the moment the Hamiltonian for the transverse phonons, 
^piion' to study the validity of the quasi ID approxima- 
tion. It is given by 



rr± ^ Eld. 



Zm ' J J 



3m: 



167reo ^ 



1^0 _ ^0|5 



(47) 



Note that for fluctuations in z-direction we have included 
the anisotropy of the dipole-dipole interactions and ob- 
tained and additional factor of 3 (see second line of 
Eq. (|T7)) ) which lifts the degeneracy of the two transverse 
phonon branches. 

The leading contribution of iJjnt is linear in the posi- 
tion operators Xi and is given by 



^1 

87reo 4" 



r0|5 



{Xi - Xj)® 



(48) 



Transverse phonons couple to excitons only in second or- 
der of yi and Zi and can be omitted in a first approxima- 
tion. 

Our goal in the following is, to study the properties of 
excitons and phonons in a MDC in the presence of an 
additional harmonic confinement potential and compare 
it with the homogeneous crystal studied in Sec. IIIII Since 
for an inhomogeneous system momentum is no longer 
a good quantum number we introduce a general set of 
exciton eigenstates, 



|l„)^i?i|Oe)=^C„(*)| 



.91 ■ 



■9n) 



(49) 



where for n = 1 . . . the C„(i) are the normalized coef- 
ficients of the eigenstates of 7?oxc and the correspond- 
ing set of (approximate) bosonic creation operators. In 
analogy, we express longitudinal displacement operators 
Xi in terms of phonon annihilation (creation) operators 
an {ai), 




^ \j 2mu>(n) 



(50) 
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with Cn{i) the mode function for phonons of frequency 
u(n) . Ignoring transverse phonons for the moment, the 
total Hamiltonian for the inhomogeneous dipolar crystal 
is then of the form 



i^MDC = J2 E{n)RiRn + 

n rn 

+ M{m,n,n'){ara + al^)RlRn' . 



(51) 




and we are left with the evaluation of the exciton spec- 
trum E{n), the phonon spectrum uj(m) and the coupling 
matrix elements M(m,n,n'). Although in the inhomo- 
geneous case we cannot derive exact analytic expressions 
for those quantities, good approximate solutions can be 
found in the limit of large N. The derivations of the fol- 
lowing results which are summarized in App. [B] an d [Cl 
are based on a similar calculation by Morigi et al. [33| 
for the phonon spectrum of a harmonically confined ion 
crystal. 



1. Density profile 

In a first step we need to determine the equilibrium 
positions of N dipoles confined by the harmonic trapping 
potential, Vt{x) = To do so we define the 

density n(x^) = l/jx^^^^^i*! which in the limit of large N 
and a slow variation of the trapping potential becomes a 
continuous function n{x). In this limit the total potential 
energy of the crystal is 



pot 



dx 



2 2 

mv X 



5An{x)+a-i)Csn^{x), (52) 



with C3 = /ig/(47reo) and 5c the chemical potential to fix 



the particle number, 
the density 



The energy E'pot is minimized for 



n{x) = 71(0)^1-4x2/^2^ 



(53) 



where the density at the center of the trap n(0) and the 
length of the crystal, L, are given by 



n(0) 



^2/5 



2C(3)i/^ 
L = A X N X ao , 




(54) 
(55) 



Here we introduced the numerical constant A = 
5r(5/6)/(r(l/3)V^) ~ 1.19. In Fig. [5] we compare the 
analytic result for n{x) with a numerical evaluation of 
the equilibrium positions x^ for molecule numbers up to 
N = 1000. We find excellent agreement between nu- 
merical and analytic results even for a small number of 
molecules. 

Eqs. ([M]) and ([55)1 express the dependence of the cen- 
ter density and the length of the crystal as a function of 
the molecule number N and a given trapping frequency 



Figure 6: a) Density profile n{x) for a harmonically con- 
fined dipolar crystal for A'^ = 100 plotted in units of a = 
{mv'^ /CaY^^ and n = 1/a. b) Length L [solid line) and cen- 
ter density n(0) [dashed line) are plotted as a function of the 
particle number TV. In both figures analytic solutions [lines) 
derived in Sec. IV B II are compared a numerical calculations 
[crosses). 



V. For comparison with the case of a homogeneous crys- 
tal discussed in Sec. IIIII it is more convenient to fix the 
density at the trap center, n(0), and adjust the trap- 
ping frequency v accordingly. In the following we express 
quantities in units of ag = l/n(0) and the corresponding 
energy scale Udd = Mg/(47''eoao). The dimensionless pa- 
rameter 7 = Udd/{fi'^ /maQ) then gives the ratio between 
potential and kinetic energy at the center of the trap. 
Note that for a fixed number of molecules we can use the 
relation 



2-5/2 ^ y^AVC(3) xhiy^ Udd/N , 



(56) 



to switch between the energy scales of the trapping fre- 
quency and the dipole-dipole interaction. 



2. Exciton spectrum 

Based on set of equilibrium positions a;^ given by the 
density profile n[x) derived above we now evaluate the 
energy spectrum of the exciton Hamiltonian i/cxc (|45p . 
As we show in App. [B] in the long wavelength limit and 
omitting the energy offset huj^g the exciton spectrum has 
the form 



E{n) 
nUdd 



2C(3)- Aj5jv+log 



N 



1/2 



n - 1/2 
N 



with numerical constants A - 
3 + log(A^log(iV/2)/32C(3)). 



(57) 

4yC(3)/A and Bn = 
Similar as for a homo- 
geneous crystal the long range character of the dipole- 
dipole interactions leads to logarithmic corrections com- 
pared to a harmonic spectrum which would result from 
nearest neighbor interactions. The corresponding exci- 
ton modefunctions C„(a;) are in a good approximation 
given by C„(a;) ~ <^n(x,an) with $„(a;, cr„) the standard 
harmonic oscillator eigenfuctions defined in Eq. (|B9p and 
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is of the form 



Figure 7; Exciton spectrum E{n) for a harmonically confined 
MDC for A'^ = 800 molecules. Numerical results {solid line) 
are compared with analytic approximations derived in App.lBl 
for the long wavelength (LW) and short wavelength (SW) 
limit {dashed lines). The dotted line indicates the LW result 
for the simply harmonic approximation given in Eq. (|B10|I in 
App.E 



an a mode dependent width 



7V- 
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1/2 



(58) 



In the short wavelength hmit the exciton energies de- 
crease linearly with increasing n down to a minimum 
energy Smin = — 3C(3)/2 x KUdd- The correspond- 
ing, rapidly oscillating modefunctions can be written as 
C„(a;0) = {-l)Cnix°) with (7„(a;°) ~ $^_„+i(a:,a) a 
slowly varying envelop function of width a ~ 0.61 \/N. 
Since the fast oscillations in the short wave length limit 
wash out the effect of long range interactions the spec- 
trum is purely harmonic. The analytic results are com- 
pared with a numerical diagonalization of i?cxc in Fig- El 
Note that the total the width of the exciton spectrum, 
AE = 7({3)KUdd/4: is exactly the same as in the homo- 
geneous case. This is due to the fact that both, long 
and short wavelength excitons are located at the center 
of trap where the density is almost constant. 



Lu{m) ~ lyy/l + (Sto^ — m — 2) /2 . 



(59) 



We recover the exact results for the center of mass mode 
uj{l) = V and the breathing mode a;(2) = \/hv and obtain 
a roughly linear phonon spectrum a;(m) ~ x 1.22 x m 
for larger m. The corresponding modefunctions Cm(i) 
describe collective oscillations extended over the whole 
crystal. 

In the short wavelength limit the spectrum of longitu- 
dinal phonos is given by 



Lo{m) — lud 



40C(3) (m-l-1/2) 



31C(5)A2 



N 



(60) 



with fh — N ~ m and a Debye frequency lod = vNK x 
^93C(5)/64^(3). When we reexpress the Debye fre- 
quency in units of Udd using Eq. (|56p we find that 



ujd exactly matches the Debye frequency of a homoge- 
neous crystal given in Eq. (|24p . This apparent coinci- 
dent is based on the fact that in the short wavelength 
limit phonon modefunctions are of the form Cm(i) ~ 
(— l)*Cm(i) with an envelop function Cfa{i) that is local- 
ized at the center of the trap (see Eq. (|C9p in App. [C| is 
therefore not sensitive to the variation of the density at 
the edges of the crystal. 

In summary we find that the longitudinal phonon spec- 
trum of a harmonically confined MDC is to good approx- 
imation linear and it is hardly affected by the long range 
character of dipole-dipole interactions. The numerical 
results plotted in Fig. [5] agree well with our analytic ex- 
pressions given in Eq. ([55)1 and Eq. (pO]) and show that 
for most purposes we can simply approximate the spec- 
trum by Lo{m) ~ X m/N where the Debye frequency 
lud is the same as in a homogeneous crystal of density 
n = n{0). 



3. Phonon spectrum 

The phonon spectrum of a quasi ID MDC with a 
strong transverse trapping frequency u± and a weak lon- 
gitudinal trapping frequency v consists of one longitu- 
dinal ('acoustic') and two transverse ('optical') phonon 
branches which are plotted in Fig. [51 In App. [U] we de- 
rive approximate analytic expressions for the long and 
short wavelength limits of the individual branches. For 
the moment we restrict our discussion to the longitudinal 
modes and come back to some properties of transverse 
phonons in Sec. IV CI 

In the long wavelength limit the spectrum of longitu- 
dinal phonons determined by hK^^^ defined in Eq. (I46p 



4- Exciton-phonon interaction 

The exciton-phonon interaction Hamiltonian iJint de- 
fined in Eq. (|48l) describes transitions between excitons 
in mode C„(i) and C„'(i) by simultaneously emitting or 
absorbing phonons in mode Cm(i). Due to the absence of 
momentum conservation for a dipolar crystal in a trap, 
the transition matrix element M{m, n, n') for this process 
depends on all three indices and has a more complicated 
structure as in the homogeneous case. Using the approx- 
imate phonon spectrum a; (to) lud m/N with cj/^ given 
in Eq. ([24]) we can write transition matrix elements as 

M{m, n, n') = ~ x M{m, n, n') , (61) 
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Figure 8: Phonon spectrum of a harmonically confined MDC 
under quasi ID trapping conditions for A'^ = 800 molecules. 
The solid lines show the numerically evaluated spectrum for 
longitudinal (w) and transverse {'-iJy,z) phonons in units of Udd 
and for 7 = 30 and hi/ ± /Udd = 1-2. The dashed lines show 
the analytic results derived in App.[C]for the long wavelength 
(LW) and short wavelength (SW) limits of the 3 different 
branches. 



where in terms of normalized equilibrium positions = 
/oq the dimensionless matrix element is 

J 27 n^lyr^ x°~x° 

Mim,n,n = n -\/ — -> r-n tTttX 

^ ' ' ^ y 62C(5) V m 2 ^ - (52) 

X (c™(i) - Cynii)) (C„'(*)Cn(j) + Cn'{i)Cn{i)) ■ 

By replacing the coefficients C„(i) (cm(i)) by the contin- 
uous modefunctions C„(x) (c,„(a;)) derived in App. |B] 
and [C] it is possible to study some general proper- 
ties for the matrix elements M(m,n, n'). For exam- 
ple, for long wavelength phonons (m/N <C 1) we find 
Cm(*) — Cm{j) ~ c'„^{x^) ^ m/N and recover the scaling 
M{m,n,n') ~ y^m/N in analogy to the homogeneous 
crystal. We here do not go further into analytic details 
of M(m, n, n') and instead we use in our calculations be- 
low numerically evaluated values for A4{m, n, n'). 



C. Stability of the crystalline phase 

Our analysis so far has been based on the assumption 
that the molecules form a linear crystal with small fluctu- 
ations around equilibrium positions. This assumption is 
valid in the limit of 7 3> 1, low temperatures, kjjT <ti Udd, 
and strong transverse confinement S> Udd- Since in 
a real experiment none of these conditions is strictly sat- 
isfied we now study the stability of our system for finite 
values of 7, T, v±. We identify three processes which 
destabilize the dipolar crystal. First, longitudinal fluc- 
tuations of the molecules eventually lead to a melting of 
the crystalline structure. For T = and a homogeneous 
system this crossover has been studied numerically in 



Ref . [19|, |20| , but we are not aware of similar studies for fi- 
nite T or for finite trapping potential. Second, for a weak 
transverse confinement there is a regime where the linear 
chain is no longer the correct ground state and molecules 
order in a zig-zag configuration. This so-called 'zig-zag 
instability' is well known for a linear chain of trap ped 
ions where it has been analyzed theoretically [s^, [34. l35l| 
and verified experimentallv [36 ) ] . A third process which 
has been studied in Ref. jlTI. Il8|| is quantum tunnelling 
of molecules into regions where dipole-dipole interactions 
are attractive. 



1. Longitudinal stability of a dipolar crystal 

We first study the melting of a dipolar crystal due 
to longitudinal quantum and thermal fluctuations of the 
molecules. A rigorous treatment of this problem would 
in principle require to take into account the full quan- 
tum many body theory for our ID dipolar system which 
is beyond the scope of the present work. Instead we here 
present a much simpler calculation assuming the valid- 
ity of the phonon Hamiltonian i?pi,o,j and determine the 
parameter regime where our model is self-consistent, i.e. 
where fluctuations are small compared to the mean sep- 
aration of the molecules. 

To study local fluctuations of the molecules in the pres- 
ence of a longitudinal trapping potential we introduce the 
position dependent 'Lindcmann parameter' 

TL{x,T)^n{x)^Ax,T), (63) 

with (Aa;)^(x°,T) = {{xi+i — Xi)"^) and the average is 
taken with respect to the phonon equilibrium density op- 
erator at temperature T. By employing the sound wave 
approximation u}(m) = ujd x m/N we can write Tl{x, T) 
as 

with a universal function F[(^,t) deflned in App. [Dl in 
Eq. (jD4[) . In Fig.[9]we plot the dependence of F{£^, t) on 
temperature and position. For the homogeneous crystal 
we obtain an analogous result with F(^, t) in Eq. ((64|) 
replaced by Fhir). Due to the nearly constant density at 
the center of the trap we identify F}i{t) = F{^ = 0, r). 
From Eq. ((64|) we flnd that fluctuations are only weakly 
suppressed with increasing 7 which is already reflected in 
the 7"^/** dependence of the exciton- phonon interaction 
(see Eq. dm). 

The minimal criterion for the local stability of the crys- 
talline phase and therefore the self-consistency of our 
model is Tl{x,T) < 1. To improve this criterion we 
argue as follows. By Quantum Monte Carlo simulations 
it has been predicted [3, [l^l that at T = a crystalline 
phase appears for 7 > 7c « 1, or in terms of our Lin- 
dcmann parameter rL(0,0) < ^'(0,0)77^" ~ 0.42. It 
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Figure 9: Numerical results for the function r) intro- 
duced in Eq. ((M]) for iV = 800. a) F(^,r) is plotted as a 
function of ^ = 2x/L for different values of r = ^ kbT /Udd- 
b) Fh(r) = — 0, r) is plotted as a function of r (solid line). 
The two dashed lines indicate the two limits Fhir ^ 1) ~ 0.42 
and Fh{T 2> 1) 0.28^^. See text for more details. 



is reasonable to assume that a generalization of this cri- 
terion, rL{x,T) < 0.42, should also provide a good es- 
timate for the local existence of a crystalline phase for 
a finite temperature and for an inhomogeneous density 
profile. Note that for the parameter regime 7 « 10 — 100 
this criterion translates into F{^,t) < 0.42 x 7^/** « 1. 
From Fig. [9] a) we find that for r = fluctuations in the 
crystal are roughly constant and almost the whole crys- 
tal is in a crystalline phase. For r > fluctuations at 
the edges of the chain quickly start to increase and at 
around t 5 — 10 already a significant fraction of the 
system does not fulfill stability criterion. Therefore for 
the inhomogeneous crystal we deduce the temperature 
limit ksT < hUddj \/^- For the homogenous crystal we 
find Fhir » 1) ~ 0.28v^ (see App. |D]) which leads to a 
similar limit for the temperature of khT < 2.4 x Udd- 



2. Transverse stability of a dipolar crystal 



closed form as 

Here ctj, = 1, = 3 and f{q) is given in Eq. (|A8[) . From 
Eq. (pS)) we see that transversal phonons are 'optical' 
phonons with a offset frequency v± and in the quasi ID 
limit, hv ^ Udd they from a flat band, uiy .^{q) « iy±. 
However, when hi^± and the dipole-dipole interaction Udd 
are of the same order we find a significant reduction of 
^x.yil) foi' short wavelength phonons (see Fig. [5]). In 
particular the transverse phonons along z-direction reach 
a minimum frequency of 

_L r 279C(5) fUddY 

The meaning of ui^^^ — > is that the linear chain of 
molecules becomes a metastable configuration and a new 
ground state appears where molecules are ordered in a 
'zig-zag' configuration. To avoid this 'zig-zag' instability 
the minimum requirement for the transverse confinement 
frequency is 

Hiy^ > 6.08 X Udd/ Vl- (67) 

While the 'zig-zag' instability is a purely classical ef- 
fect there is also a quantum mechanical instability of a 
quasi ID dipolar crystal due to tunnelling events into the 
region of attractive dipole-dipole interactions. For a de- 
tailed description of this process the reader is referred to 
Ref. [l^l where this tunnelling rate has been analyzed for 
the case of two molecules approaching each other under 
ID or 2D trapping conditions. When we adopt the re- 
sult derived in Ref. [ll|. Sec. Ill A 4 and use as the 
attempt frequency we a obtain tunnelling rate 

Ftun 0^ LOU exp (-c^fi^wjWM)) , (68) 

with a numerical constant c sa 5.8. When we assume that 
the inequality (|67p is satisfied we find that the tunnelling 
rate is bound from above by Ftun ^ exp(— c^/7). We 
conclude that when the system is in a crystalline regime, 
7 > 1, and condition ([57|) is fulfilled the crystal is also 
stable with respect to tunnelling events. 



Position fluctuations of the molecules in transverse di- 
rections are described by the Hamiltonian H:^Yion given 
in Eq. (|47p . The resulting phonon spectrum of the two 
transverse phonon branches is plotted in Fig. [8] and ap- 
proximate analytic expressions for the long and short 
wavelength limit can be found in App [C] in Eqs. ICllI 
and lC12l As transverse phonons are hardly influenced by 
the longitudinal trapping potential we use in the follow- 
ing discussion the exact spectrum of transverse phonons 
in a homogeneous crystal, which can be written in a 



3. Summary 

In summary we find that apart from 7 > 1 the stability 
of a homogeneous quasi ID dipolar crystal requires the 
following relation between the relevant energy scales, 



max 



ma, 



,0.42 keT 







<Udd<^/l 



6l)8 



(69) 



In the inhomogeneous case the bound on the temperature 
should be replaced by fc^T < 5 Udd/ \/l to guarantee a 
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crystalline structure over a large fraction of the system. 
Finally we emphasis that numerical values for the bound 
on the temperature should not be considered as precise 
numbers. One hand this is due to our simplified model 
we used to calculate these numbers. On the other hand in 
a ID configuration we anyway expect a smooth crossover 
from a crystalline to a liquid regime with no precisely 
defined transition point. In that sense Eq. (|69p deter- 
mines the parameter regime where our model Hamilto- 
nian Hmdc (EH) is valid. Outside this regime we do not 
expect a complete break down of our model but higher 
order corrections should be taken into account. 



D. State transfer fidelities for spin ensemble qubits 

So far in Sec. |V] we have discussed several aspects of a 
MDC in the presence of a finite longitudinal and trans- 
verse confinement potential and we have identified the 
stability criteria for such a system. In particular we 
have focused on a specific scenario where molecules are 
trapped in electrostatic potentials which provides a con- 
venient way to achieve strong transverse confinement and 
the possibility to bring molecules close the surface of the 
cavity electrode to enhance the coupling strength. How- 
ever, the severe restriction on trappable rotational states 
limits the induced dipole moments to fig ^ 0.2/io ~ ID 
while at the same time resonant dipole-dipole interac- 
tion and therefore exciton-phonon interactions are quite 
large (k ~ 10). In this section we finally want to show 
that even under those unfavorable conditions we may still 
achieve high fidelity state transfer operations between the 
cavity and ensemble qubits encoded in collective spin ex- 
citations. 



1. Ensemble-cavity coupling 

In order to couple the two spin states \g) and |s) defined 
in Sec. lV Al we consider a Raman type setup where \g) and 
the rotational state |e) are coupled to the cavity field and 
|e) is in turn coupled to state |s) by a classical microwave 
field of orthogonal polarization. For this configuration 
the cavity-molecule interaction is 

ifcav-moi(i) = ^g(a;)|g,)(e,|ct + ^|s,)(e,| + H.c. , 

(70) 

with n{t) the controllable Rabi frequency of the exter- 
nal microwave field. In Eq. (j70p we have generalize our 
model to a non-uniform cavity modefunction by intro- 
ducing the position dependent single molecule coupling 
g{x) — gu{x). Here g is the maximum coupling constant 
at the center of the trap and < u(a::) < 1. Note that 
u{x) varies at least on the scale of the cavity wavelength 
Ac, but may in principle be designed to have an arbitrary 
shape. Under the two-photon resonance condition and a 
detuning A = ojc — '^eg ^ ff's/^i Q,{t) we can eliminate 



the excited state and obtain 

Jfcav-mol(t) - gR{t) (^jc + S,) , (71) 

with a Raman coupling strength gR{t) — gNCl{t) /2A. 
Note that due to the non-uniform cavity-molecule cou- 
pling we obtain a modified ensemble coupling g^ = 
gV^cS with an effective number of molecules N^h — 
J \u{x)\'^n{x)dx. Similarly, ensemble qubit states are de- 
fined as 

|1.) ^ 5l|0,) = Yl ■ ■ ■ • . . 5n) • (72) 

Assuming for simplicity f2(i)/2 « g^ the Raman cou- 
pling Hamiltonian (|7T|) provides the basic ingredient for 
a swap operation between the state of the cavity and the 
spin ensemble qubit in a time Tg — TTA/2gfj ^ ^/gN- For 
a realistic estimate of g^ we consider the predictions for 
the single molecule coupling strength g which are given 
in Ref. 1] for CaBr with = 4.3 D. With d the distance 
between the molecules and the cavity electrode we obtain 

gAr/(27r) ~ 40kHz x y/N^/d[fj.m] . (73) 

As in the crystalline phase we do not have to care about 
motional diffusion of molecules the length of the crystal 
can in principle be as large as L w Ac/ 2 meaning that 
even for a ID crystal the number of molecules can be as 
high as Ncff ~ Ac/(2Aao) ~ 10^. Using a moderate trap 
distance d — 0.5 /im we end up with a collective coupling 
strength in the order of f^^v / (27r) w 25 MHz, which can in 
principle be pushed into the 100 MHz regime by going to 
trap-surface distances of ci ~ 0.1 /xm. 

2. Decay of rotational ensemble qubits 

As decoherence processes for spin ensemble qubits dur- 
ing gate operations are due to a finite population of the 
rotationally excited state |e) we study in a first step the 
decay of a rotational ensemble qubit state |le). Follow- 
ing the calculations of Sec. IIII Gl the initial decay of the 
excited state probability is Pe{t) ~ 1 — where com- 

pared to the homogeneous case the decay rate W now 
contains two contributions, 

W'=(^^y (1,^, + ^^ . (74) 

The first term in Eq. ([74|) is the exciton dispersion 
which arises in an inhomogeneous system from the mis- 
match between cavity modefunction and exciton eigen- 
functions. In terms of the normalized exciton spec- 
trum E{n) = {E(n) — huj eg) / {nU dd) and the overlap 
Zn = Y.i C^n{i)u{x'-) it is given by 

Icxc = ^^'(7^)4 - {Y.Ein)ziy. (75) 

n n 
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For a cavity modefunction u{x) = cos(7rx/Ac) and L < 
Ac/2 we obtain numerical values of Xoxc — 0.40 — 0.11, 
where the lower value corresponds to L — Ac/2. Of 
course, a more sophisticated trap or cavity design would 
reduce this value even further towards Xoxc — 0, e.g., in 
the limit of a flat bottom trap. 

The second term in Eq. ([7^ describes the decay due 
to exciton-phonon interactions and is defined as 

2'phon= Yj ZnZn'\M{m,n,n')f {2N{Lu{m)) + l) , 

(76) 

with normalized coupling matrix elements A4(rn,n,n') 
defined in Eq. (|62p . A numerical evaluation of 2phon for 
the case L <C Ac shows that in the zero temperature 
limit Iphon ~ 1-38 while for high temperatures we obtain 
2phon ~ 11-3 X T, with T = kBT/Udd- The crossover 
point between the two regimes is r « 1. 



3. State transfer fidelities for spin ensemble qubits 

For a simplified discussion of the state transfer fidelity 
between the microwave cavity and spin ensemble qubits 
we identify the gate fidelity J- with the probability to 
convert a single cavity photon, |lc), into a spin excitation 
|ls). The fidelity is degraded by two processes. First, 
the spin state jig) decays due to the finite admixture 
of the rotational ensemble state jig). For a detuning 
A » l-ffMDcl and the fast gate times Tq considered in 
this paper this decay is quadratic and the corresponding 
rate Ws ~ (gN/^)"^ x is proportional to W defined 
in Eq. ([74|l . but suppressed by {gM/^Y- Second, the 
photon state |lc) decays linearly with the cavity decay 
rate Fc. Assuming that each of the two processes acts 
for approximately half of the gate time Tq we obtain a 
total gate fidelity 



1 



ttVF 
4A 



ttFcA 
1^ 



(77) 



We optimize the gate fidelity for a detuning : 
(Trg^VF^/Fc)^/^ which results in a maximal fidelity of 



3 /TT 

1-4(2 



TcW 
9% 



(78) 



Discussion. We now consider a specific example us- 
ing the molecule CaBr with /io — 4.3 D, to = 120 amu 
and B/27r = 2.8 GHz. From Table |T] a) we find that at 
the sweet spot the induced dipole moment is fig ~ 0.7 
D and k — 10.5. To achieve a stable crystal we choose 
a lattice spacing of aq — 70 nm which corresponds to 
7 « 13 and C/dti/(27r) = 215 kHz. From our estimates 
on the stability of the crystal in Sec. IV CI we obtain the 
conditions fc^T < 2Udd « 20 /xK and i^±/2tt > 360 kHz. 
Both requirements seem feasible with on-chip cooling and 
trapping techniques proposed in Ref. [ij. Trapping a 



moderate number of molecules N « 10^ at a distance 
d 0.5 fim above the cavity electrode we obtain a col- 
lective coupling strength of gN/'^'^ ~ 8 MHz. Using a 
superconducting microwave cavity with a quality factor 
Q « 10^ as demonstrated in Ref. 15] the decay rate of 
the cavity is as low as Fc — ojc/Q ~ 27r x 10 kHz. With- 
out any mode matching but assuming temperatures of 
ksT < Udd/^/l w 3 /xK we obtain W/{2n) « 2 MHz. 
Inserting these values into Eq. ((78| we obtain a gate fi- 
dehty of F ~ 0.994 and a gate time of Tq ^ 0.14 ^s. For 
a more optimistic choice of parameters with gpf/2'K — 25 
MHz and oq = 100 nm (which corresponds to 7 sa 9, 
Udd/'^T^ ~ 75 kHz and a temperature requirement fc^T < 
1 /iK) we immediately obtain gate errors of below 10"'^ 
at even shorter gate times. 

In conclusion we find that while in an electrostatic 
trapping configuration exciton-phonon interactions are 
quite high, it can be overcome by the high collective cou- 
pling strength due to the high densities in the crystalline 
phase and the low trap surface distance. However, for 
further improvements it might be necessary to consider 
magnetic trapping techniques where exciton-phonon in- 
teractions can be highly reduced while at the same time 
temperature requirements would be less stringent. A 
second interesting alternative is the choice of rotational 
states \g) and |e) which are given in Table HI example 
e). Both states are weak field seekers and at the same 
time they satisfy the decoupling condition k + e = 0. 
However, as molecules in state \g) and |e) have a differ- 
ent induced dipole moment and feel a different trapping 
potential, this configuration would require a flat bottom 
trap in longitudinal direction and ground state cooling in 
transverse directions. 



VI. SUMMARY & CONCLUSION 

In this paper we have investigated the storage of quan- 
tum information encoded in collective excitations (en- 
semble qubits) of long-lived rotational or spin degrees 
of freedom in a self-assembled dipolar crystal of polar 
molecules. This provides a high fidelity quantum mem- 
ory which can be coupled to a superconducting strip line 
cavity which in the spirit of Cavity QED provides a cou- 
pling to a solid state quantum processor. The main re- 
sults are summarized as follows. 

In the first part of this work we have studied the dy- 
namics of rotational excitations (= excitons) in a self- 
assembled molecular dipolar crystal (MDC) which maps 
to a polaron type model with excitons interacting with 
the phonon modes of the crystal. While in general the 
exciton-phonon interactions plays the dominant role as a 
decoherence mechanism in this system, leading to a decay 
of rotational ensemble qubits, we have identified certain 
'magic' configurations where long wavelength excitons - 
which includes the ensemble qubit state - decouple from 
the phonon modes. Furthermore, quantum information 
encoded in spin ensemble qubits is naturally protected 
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from dipole-dipole interactions, and the exciton-phonon 
interactions affects the ensemble quantum memory only 
during gate operations. 

In the second part of this paper we have studied in 
detail a specific scenario with molecules trapped in elec- 
trostatic potentials (e.g. an on-chip electric trap) and 
quantum information encoded in collective spin excita- 
tions. We have discussed modification of the exciton and 
phonon spectrum due to the presence of a longitudinal 
trapping potential and analyzed the stability of a MDC in 
a quasi- ID geometry. An estimate of the expected state 
transfer fidelities between the microwave cavity and spin 
ensemble qubit for this setup shows that under reason- 
able experimental conditions fidelities of > 0.99 can be 
achieved for a total gate time well below 1 jis. Optimized 
conditions would result in gate errors of ~ 10"'* which 
would allow fault tolerant quantum computing (37| . This 
specific example demonstrates the potential of MDCs in 
the context of hybrid quantum computing since high gate 
fidelities and long storage times are combined with gate 
times that are compatible with decoherence times scales 
in solid state based quantum computing. 
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Appendix A: HOMOGENEOUS DIPOLAR 
CRYSTAL 

In this Appendix we briefly summarize the deriva- 
tion of the exciton spectrum i?(k), the phonon spec- 
trum LU\ (k) and the coupling matrix elements Mx (q, k) 
for a homogeneous dipolar crystal in ID and 2D. All re- 
sults in this appendix are expressed in units of the lat- 
tice spacing oq and the corresponding dipole-dipole en- 
ergy Udd = A*g/(47reoao)- In these units the parameter 
7 = Udd/{^ /maQ) plays the role of the dimensionless 
mass of the molecules. 

Excitons. We start with the exciton Hamiltonian -ffexc 
as defined in Eq. ([T5]). In dimensionless units it is given 
by 



9 51 Ir-O 



2^ iry 



r-0|3 



(Al) 



with Eeg = hiOeg/Udd and Kij defined in Eq. ((T3)) . In 
the limit of a low number of rotational excitations we 
can expressing iJcxc in terms of the exciton operators 
R]^ — e*'"'i \ei){gi\ by making the substitutions 



N 

k,k' 

+ K{e+^'^°e-'^''" +e+^''"e'■^^''' 



and — X^k^k^k- Evaluating the resulting 

expressions we end up with a diagonal Hamiltonian of the 
form iJcxc = X^k ^(l^)^k^kj where the energy spectrum 
i?(k) = Ecg + eJiO) + KJ(k) is given in terms of the 
dimensionless function 



J(k) = 



cos(kr°) 



r0|3 



(A3) 



For the ID crystal with ri = i we can evaluate J{k) and 
obtain 

j(fc)^2fj^^ =Li3(e-^-)+Li3(e^'^), (A4) 

with Li„(z) the polylogarithm function. In 2D the equi- 
librium positions form a triangular lattice with basic 
lattice vectors ai = (1,0) and 02 = (l,\/3)/2 and we 
evaluate the function J(k) numerically (see Fig. |4]). By 
replacing the summation in Eq. (jA3[) by an integral we 
obtain the linear behavior, J(|k| — > 0) — J(0) ~ |k|, for 
the long wavelength limit. For a more detailed study 
of the spectrum of (rotational) excitons in 2D analytical 
tools developed in the field of 2D Wigner crystals [3g| or 
2D spin waves ||39| can be applied to handle the slowly 
convergent sums in Eq. (jA3p . This analysis will be the 
subject of future work. 

Phonons. In a next step we consider the Hamiltonian 
of longitudinal phonons, iJphon, which is given by 



phoii 



E 
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3^ 5 [(x,, - Xj) ■ n°.]^ - [x, - Xj]^ 

4 |j.0_j.0|5 



(A5) 

As -ffphon is quadratic in position and momentum opera- 
tors we can rewrite it terms of phonon annihilation and 
creation operators, i/phon = I]q,A ^^A(q)a^(q)aA(q). To 
find the phonon spectrum a;A(q) we change into the 
Heisenberg picture and make the ansatz 
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where in the 2D case the vectors eA(q) are the two 
orthonormal polarization vectors of the two phonon 
branches. This ansatz in combination with the Heisen- 
berg equations Xi{t) = — [iJphon, Pi]/l leads to eigenvalue 
equation 



'^A(q)eA(q) =y^(q)eA(q) 



(A7) 



Here -4(q) is a single valued function in ID and a 2 x 2 
matrix in 2D. For the ID case we rewrite it as A{q) — 
f\q)h with 
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Sin2(qj72) 



(A8) 



= 12[2C(5)-Li5(e''^)-Li5(e-^')] 
For the 2D crystal the matrix A{c\) is defined as 



^(q) 
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(1 



^015 



(A9) 

By numerically solving the eigenvalue problem for the 
matrix ^(q) for each value of q we obtain the two phonon 
branches liJ\{<\) and the corresponding polarization vec- 
tors e\. In analogy to the ID case we express the re- 
sulting phonon spectrum in terms of the two rescaled 
functions /A(q) = ujx{q)^ (see Fig.S]). 

Exciton-phonon interactions. Finally, we consider the 
first order exciton-phonon interaction, -ffint, given by 



int 



r0|5 



(AlO) 



Using Eqs. (|A2p and (|A6p we reexpress the interaction 
Hamiltonian in terms of exciton and phonon operators 
we leads us to the form of Hint given in Eq. (fTOj) with a 
coupling matrix element 



MA(q,k) 
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r^^eA(q) / 
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By introducing the function 

3 ^r"-eA(q) 

ffA(q) = ^E^ir 



sm 



(qr°) 



(A12) 



we can bring Mx (q, k) into the form given in Eq. ([25 
For ID we can evaluate the sum analytically. 
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(A13) 

In 2D the evaluation of the functions 5a (q) requires a 
numerical summation of the right hand side of Eq. (|A12[) . 



Appendix B: EXCITON SPECTRUM OF A 
HARMONICALLY CONFINED MDC 



In this appendix we derive approximate analytic ex- 
pressions for the exciton spectrum in a harmonically con- 
fined dipolar crystal in ID. All quantities in this appendix 
are expressed in units of the mean molecule separation 
at the center of the trap a(0) = l/n(0) and the corre- 
sponding dipole-dipole energy Udd (see Sec. IVB Tl) . 

We start with the exciton Hamiltonian i?exc given 
in Sec. IV Bl in Eq. ([45]) and look at the eigenvalue 
equation E{n)R\ — [TJexc-Rjj] for exciton operators 
Ri = J2^ Cn{i)\ei){gi\ with n = 1...N. In the limit of 
low number rotational excitions the resulting eigenvalue 
equation is of the form 



C'njj) 
_ ™0|3 



(Bl) 



with Eeg ~ huJeg/Udd- For convenience we omit this con- 
stant energy offset in the following calculations. Our goal 
is to derive approximate analytic solutions of Eq. (|Bip in 
the long and short wavelength limit. 

Long wavelength limit. For a large number of molecules 
N the density n(x) varies slowly over the extension of the 
crystal and we make the approximation 



E' 
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r0|3 



2C(3)n3(a:°) , 



(B2) 



to convert Eq. (|Bip into the form 



3^1 '""J 



|X0-X0|3 



(B3) 

In the long wavelength limit C„(a;^) = C„(i) can be 
approximated by a slowly varying continuous function 
C„(x) and 

C„(j) - C„(f) ^ C'^{x%x] - x\) ^ C;:(x«)(x° - x^)V2 . 

(B4) 
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Inserting this expansion into the right hand side of 
Eq. (|B3p we obtain two contributions. The first one 
proportional to the first derivative C^{x) vanishes due 
to the summation over an equal number of terms with 
positive and negative sign. In contrast, the second con- 
tribution proportional C"(x) diverges as ~ log(iV). This 
divergence is a consequence of the long range character 
of dipole-dipole interactions and means that Eq. (|B3p 
is sensitive to the spatial extension of the wavefunction 
C„(x) and the expansion in Eq. (jB4p will fail to predict 
accurate results. To handle this difficulty we proceed as 
follows. In a first step we approximate the divergent sum 

by 



(B5) 



where in a zeroth order approximation we set Eq — 
log(iV/2). Under this approximation we find that (see be- 
low) wavefunctions C„(a;) are harmonic oscillator eigen- 
functions. In a second step we use the spatial dependence 
of these zeroth order eigenfunctions for a more accurate 
reevaluation of the right hand side of Eq. (|B3p . 

Using the ansatz in Eq. (jB5P we can convert the eigen- 
value equation (IB3|) into the differential equation, 

[E{n) ~ 1C,{?>)Kr?(x)\ Cn{x) = KY.Qn{x)C'^(x) . (B6) 

For the density profile n{x) given in Eq. (|53p and with a 
rescaled variable y = 2x/L we obtain 



(l_y2)l/3 



(B7) 



with a„ = 1 - E{n)/{2C{?,)K) and 9 = C(3)iV(2So)- In 
a final approximation we expand Eq. (jB7[) to lowest order 
in and e„ and end up with the differential equation 



C:{y) + b\an~y^)Cn{y) = 0. 



(B8) 



Solutions of this equation are C„(a;) '-^ $„(x, cr) with 
= L\/Eo/8C(3) and <I'n(a;,(7) the well-known har- 
monic oscillator functions 



shape of the wavefunction C„(x). To take into account 
this mode-dependent cutoff we use the following trick. 
First, for each modefunction Cn{x) we replace the con- 
stant So in Eq. (jBSp by E(n), which has a different value 
for each mode n. Second, to derive an accurate expres- 
sion for Ti{n) we invert the expansion of C„(x") outlined 
in Eqs. (|B3P - (jB5p and express E(n) in terms of the full 
wavefunction Cn(x) as 



E(n) 



1 



C„(x")-C„(x^) 
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(Bll) 







To evaluate Eq. (jBlip we focus on the region x\ 
and neglect small variations of the density. For the mode 
functions C„(x) we insert the harmonic oscillator func- 
tions, Cn{x) ~ <^n{x,a), derived above. Around the cen- 
ter of the trap we can further approximate these wave- 
functions by Cn{x) ~ cos(knx) (or C„(x) ^ sin(fc„x) for 
odd n) with a wavevector fc^ = (2n — 1)/(t^. For even 
modes this approximation results in 
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which can now be evaluated in a closed form in terms 
of polylogarithm functions (see App. El Eq. (|A4p V Ex- 
panding the result to lowest order in A:„ we finally obtain 



S(n) 



3 -I- log 
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By replacing Eq with E(n) in Eq. (jBlOp we end up with 
the improved exciton spectrum 



E{n) 
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N 



log(Av/log(7V/2)/32C(3)) and A 



4\/C(3)/A. The corresponding improved modefunctions 
are given by Cn{x) ^ $„(x, cr„) with a n-dependent 
width 



$„(x, cr) iJ„_i(x/o-) e 



(B9) 



with Hn{x) Hermite polynomials. Using L = AN the 
corresponding harmonic energy spectrum is 



E{n) 



2C(3) 



\/32C(3)Eo (n-1/2) 
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At this point the energy spectrum E{n) given in 
Eq. (jBlOp still depends on an inaccurately defined con- 
stant Eq ~ log(A^/2) and its shape only poorly approxi- 
mates the exact spectrum plotted in Fig. [71 As mentioned 
above this discrepancy is due to the fact that the cutoff 
radius for the sum given in Eq. (|B5p is not determined 
by the size of the system, i.e Vc = N/2, but rather by the 
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Although our derivation was based on several rather 
crude approximations we find in comparison with nu- 
merics that our analytic solutions for E{n) and Cn{x) 
provide an accurate description of the long wavelength 
behavior of the exciton spectrum and the shape of the 
eigenfunctions. Due to the similarity of the underlying 
equations this approach should also be applicable to im- 
prove the linear phonon spectrum obtained in Ref. [33| 
for a harmonically confined ion crystal. 

Short wavelength limit. For the short wavelength limit 
of the exciton spectrum we make the ansatz C„(a;^) = 
Cn{i) = {—iyCn{x^) such that C'„(x) represents a slowly 
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varying envelop function for the rapidly oscillating mode- 
function C„ (x) . Inserting this ansatz into the eigenvalue 
equation (jBl|) we can proceed in the derivation of Cn{x) 
as explained above for the long wavelength limit. How- 
ever, in the short wavelength limit the fast oscillating 
eigenmodes cancels the effect of long range interactions 
and in contrast to the divergent sum in Eq. (|B5p the cor- 
responding term in the short wavelength limit has a well 
defined value, 



■log(2)n(a;°). 



(B16) 



Therefore, the resulting spectrum in the short wavelength 
limit is purely harmonic and has the form 



E{n) = K 



3C(3) , V24C(3)log(2) (n-1/2) 
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with n = N — n + 1. The envelop functions C'n{x) ^ 
^n{x,(7) are harmonic oscillator eigenfunctions with a 
width 



7VAv/log(2)/6C(3) . 



(B18) 



Appendix C: PHONON SPECTRUM OF A 
HARMONICALLY CONFINED MDC 



for the normal modes Cm{i)- In the long wavelength limit 
we can replace the discrete set of coefficients Cm{x^) = 
Cm(*) by a continuous function Cm(x) and make the ap- 
proximation 

c™(j) 0, c,„(z) + c'ix°){x° - x°) + c"(x°)(x« - x°r/2 . 

(C3) 

As the modefunctions of long wavelength phonons extend 
over the whole length of the crystal we include variations 
of the density n{x) given in Eq. ([55)1 . i.e., 
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Keeping terms up to second order in the derivatives (x) 
and n'{x) Eq. (|C2[) transforms into the differential equa- 
tion 

nHyyLiy) + 4^ri\y)n'{y)cUy) ~ a„,c,n{y) = , (C5) 

with y — 2x/L and a™ — 2(1 — a;^(m)/z/^)/3. Here we 
made use of the identity ^v'^L'^ = 32C(3) (see Sec. IVB 1[ 
Eq. ([5^ ). For the density n(y) ~ yl — we find 
n?'{y)n'{y) ~ —2/32/ ^-nd the differential equation (jCSP 
can be solved by the ansatz Cm{y) = Xlfe '^kV^ ■ The co- 
efficients obey the recursion rule 

ak+2 — o-k [k{k - 1) + 8A:/3 + a,n] /{k + l)(fc + 2) . 



In this appendix we derive approximate analytic ex- 
pressions for the spectrum of longitudinal and transverse 
phonons in a harmonically confined MDC. Our calcula- 
tions are based on a similar approach used in Ref. (33| 
for the phonon spectrum of a ID ion crystal. However, 
in contrast to the ion crystal or the exciton spectrum de- 
rived in App.[Bl the long range character of dipole-dipole 
interactions has no severe effect on the phonon spectrum 
of a dipolar crystal which simplifies calculations. Here we 
outline only the derivation of the spectrum of longitudi- 
nal phonons. The spectrum of transverse phonons can 
be derived along the same lines and we give the results 
at the end of this appendix. Note that throughout this 
appendix we express results in units of a(0) = l/n(0) and 
Udd as defined in Sec. IVB II 

To derive the eigenspectrum uj{m) with to = 1 ... of 

the Hamiltonian ifpi^^jj given in Eq. (j46p we change into 
the Heisenberg picture where position operators Xi{t) 
obey the equation of motion. 



Xi{t) = -i^^Xi 



- - E 



Xt{t) - Xj{t) 



r0|5 



(CI) 



with D = hu /Udd the normalized trapping frequency. Us 
ing the normal mode decomposition given in Eq. ([S"" 
Eq. (jCip translates into the eigenvalue equation, 



7[^^-c.2(to)]c„,(z) = -12 5:^ 



r0|5 



(C2) 



The quantization condition for mode to, ak+2 = 0, V/j > 
TO — 1, follows from the normalizability of the resulting 
polynomials and translates into the spectrum 



w(to) = z> X + (3r 



2)/2. 



(C6) 



The corresponding modefunctions Cm{x) are polynomials 
extended over the whole length of the crystal. 

In the short wavelength limit we make the ansatz 
Cm{i) — {—^Ycm{x^) and repeat the calculations from 
above for the slowly varying envelop function Cm {x) . For 
y = 2x/L we obtain 

~„ , , _^ 4n'(y) . {f3N^n^{y) + a^) . , . ^ 



3n(y) 



n^{y) 



(C7) 



with /? = 31C(5)AV24C(3) and = 8{l-uj^{m)/C'^)/9. 
We can simplify Eq. (|C7p by neglecting the term propor- 
tional to n'{y) and by expanding the remaining equation 
up to second order in y. This approximation is valid since 
we will find below that modefunction of short wavelength 
phonons are located at the center of the chain. We end 
up with 

c'Liy) + [N^P + - {2N^P/3 ~ arn)y^] drniy) = . 

(C8) 

Solutions of this equation are of the form 

c„,{x)^H^{x/a)e'^"/^^\ (C9) 
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with TO = iV - m and = 7VAV3/(80/?). The corre- 
sponding spectrum is given by 

with a Debye frequency ujd = vN y^9(3/8. 

Transverse phonons. Along the same hues as shown 
for the longitudinal phonons we calculate analytic expres- 
sions for the spectrum of transverse phonons determined 
by the Hamiltonian H^^^^ given in Eq. (|47)l . In the long 
wavelength limit the resulting spectrum is 

± , ^ 4C(3) (3to2 - 

with ay = I, az — 3 and i>± =^ hv^/Udd- For the short 
wavelength limit we obtain a spectrum of the form 

(C12) 

with numerical constants A — 93C(5)/8 ~ 12.05 and B = 
a/5C(3)32/31C(5) ~ 2.055. 

Appendix D: LINDEMANN PARAMETER 

In this appendix we calculate the local Lindemann pa- 
rameter Tl{x, T) as defined in Eq. (|63|) both for the ho- 
mogeneous and the inhomogeneous crystal. Note that in 
the following wc express all quantities in units of oq and 
Udd- Using the normal mode decomposition of operators 



Xi given in Eq. (|A6|) for the homogeneous system and in 
Eq. ([50]) for the inhomogeneous case we obtain 

ri(x°, T)^Y1 0^^\ck{i + l)-cu{i)\\2N{u:{k)) + ^) , 

(Dl) 

with N{uj{k)) the thermal occupation number of mode k. 
For a homogenous system k is the quasi momentum, Ck{i) 
are plane waves and uj{k) — f{k)/^. In the homoge- 
neous system we replace k by the index to = 1, . . .TV 
and the normal modes Cm(i) with the corresponding 
spectrum uj{m) are discussed in Sec. IV Bl and App. [C] 
For simplicity we adopt the sound wave approximation 
uj{m) ~ Qfu/^/YX m/iV with q;d « 6.95. The Lindemann 
parameter then has the general form 

rLix,T) = ^xF{^^2x/L,T = ^T) . (D2) 
For the homogeneous crystal F(^,t) = Fii{t) with 

m 

In the two limits of interest the numerical values of this 
function are Fh{T ^ 0) ~ 0.424, and Fh{T > 1) ~ 
0.278 X y^. In the inhomogeneous case we obtain 

(D4) 

Numerical values of this function arc plotted in Fig. [51 
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